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Abstract 

We explore the Higgs-Gauge configuration space in the standard electroweak theory. 
We outline a general prescription that uses the non-trivial topology associated with 
the gauge group of the theory, to find known solutions of the Euclidean classical equa- 
tions of motion and motivate the existence of novel ones. In Minkowski spacetime 
we present evidence for the existence of approximate breathers - long-lived, spatially 
localized, temporally periodic configurations. We consider heavy fermion quantum 
fluctuations about static Higgs-Gauge conflgurations, and argue for the existence of 
stable fermionic solitons. These could resolve the fermion decoupling puzzle in chi- 
ral gauge theories. We describe our search for a fermionic soliton within a spherical 
ansatz, and discuss the quantum corrected sphaleron and the emergence of new barri- 
ers suppressing the decay of heavy fermions. Finally, we consider electroweak strings 
and how they could give rise to stable multi-quark objects. 
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Chapter 1 
Introduction 



The Standard Model of particle physics has several well-known configurations of Higgs 
and gauge fields that are solutions of the classical equations of motion. These have 
rich phenomenology associated with them. For example, the electroweak instanton 
PI |2j in Euclidean spacetime mediates non-perturbative fermion number violation 
through quantum tunneling. There are compelling reasons to expect the existence of 
novel configurations that drive physics ranging from decoupling of heavy fermions to 
electroweak baryogenesis. I explore these possibilities in this work. 

I begin by examining the space of solutions of the classical equations of motion in 
the electroweak theory, in Chap.|21 I describe my work with E. Farhi and N. Graham, 
in which we use topologically non-trivial maps into the gauge group to construct and 
motivate solutions in Euclidean spacetime. This method is a generalization of ideas 
introduced by Manton and Klinkhamer |4j. It encapsulates the various known 
solutions (electroweak instanton, sphaleron, strings, etc.) into a unified framework. 
Moreover, the use of all the topological properties of the theory leads to the possibility 
of new, unstable solutions. 

I then consider Minkowski spacetime, and argue for the existence of approximate 
breathers - long-lived, spatially localized, temporally periodic configurations - in 
the electroweak theory. E. Farhi, N. Graham and I have used the mechanism that 
creates intrinsic localized modes in anharmonic crystals [3 E] to demonstrate the 
presence of approximate breathers in the 1+1 dimensional Abelian Higgs model. We 
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are currently looking for such objects within a spherical ansatz in the electroweak 
theory. These breathers have lifetimes that are orders of magnitude larger than all 
scales in the problem and challenge the notion of naturalness in field theories. On a 
more phenomenological front, electroweak breathers could create out-of-equilibrium 
regions in space during the electroweak phase transition and drive baryogenesis. 

All known (and expected) spatially varying, static solutions in the electroweak 
theory are unstable and are generically called sphalerons (to distinguish them from 
stable solutions or solitons). They do not have any associated quantum extended- 
particle states. The discovery of a stable configuration would result in a soliton sector 
in the Hilbert space of states, in addition to the familiar vacuum sector. There is 
no topological reason for stability of a static configuration in the electroweak theory, 
but a non-topological soliton (corresponding to a local minimum of the energy) may 
still exist. However, in the absence of a topological beacon, it is difficult to search 
for such an object. If we consider quantum fluctuations around classical configura- 
tions, then there are compelling reasons to expect the existence of quantum solitons, 
and well-understood mechanisms to guide the search for them. Such objects could 
be stabilized by virtue of carrying a conserved quantum number, in analogy with 
topological solitons that carry a topological charge. 

In Chap. 121 1 explain how the quarks and leptons in the Standard Model could 
be strongly bound by certain configurations of Higgs and gauge fields, giving rise 
to the possibility of fermionic solitons. These would allow heavy fermions to de- 
couple from the theory because the lower energy fermionic solitons would carry 
their quantum numbers and maintain anomaly cancellation. Since these solitons 
are quantum-stabilized, we have to include quantum corrections to their energy when 
analyzing their stability. I review an efficient method based on scattering theory that 
allows an exact computation of the one-loop quantum corrections to the energy non- 
perturbatively, with physical on-shell renormalization carried out in the perturbative 
sector [71 IH] . This makes it feasible to carry out a variational search for fermionic 
solitons. I describe my search (with E. Farhi, N. Graham, R. L. Jaffe and H. Weigel 
jni) within a spherical ansatz and discuss the quantum-corrected sphaleron and the 
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emergence of new barriers that suppress heavy fermion decay. We find no evidence 
for a spherical fermionic sohton. However, this does not preclude the existence of 
such objects outside the ansatz. 

In Chap.m I briefiy review the family of electroweak strings PU]- These are static, 
unstable solutions of the classical equations of motion that localize energy within a 
cylindrical region of space. The Higgs condensate is suppressed in the core of a string, 
and so heavy quarks are bound along its length and resist its disintegration. If this 
mechanism could stabilize electroweak strings, then they would constitute a crucial 
ingredient in a viable scenario for electroweak baryogenesis without requiring a first- 
order phase transition [T^. Furthermore, a gas of electroweak strings would have 
negative pressure and could contribute to the dark energy required for the observed 
cosmic acceleration [T^ I13j . An analysis of the stability of these objects requires a 
computation of the fermion quantum correction to the energy. This is most efficiently 
done using scattering theory methods and there are technical challenges associated 
with the long-range nature of the potential generated by such configurations. So, as 
a stepping stone to the calculation I investigate the similar, but simpler, problem of 
one-loop quantum corrections to the energies of magnetic fiux tubes in QED (in collab- 
oration with N. Graham, M. Quandt, O. Schroeder and H. Weigel). I comment on the 
puzzles and challenges associated with the presence of fermions in such backgrounds. 
We find that when we include a region of return fiux (which is the truly physical case) 
the scattering potential becomes short-range and the calculation becomes tractable. 
As the return fiux is made more diffuse, the quantum-corrected energy of the config- 
uration has a well-defined limit and gives the energy associated with an isolated fiux 
tube. We also discover that when the energy is properly renormalized, the results 
in two and three spatial dimensions become similar. These findings suggest that a 
natural next step of research would be to examine electroweak vortices in two spatial 
dimensions by including a region of spread-out return fiux. A fermionic soliton in 
that theory would be indicative of a stable electroweak string in the Standard Model. 
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Chapter 2 



Classical Solutions 



We briefly review the bosonic sector of the electroweak theory. We use topologi- 
caUy non-trivial maps into the gauge group to explore the space of solutions to the 
classical equations of motion, in Euclidean spacetime. We show how the method 
can be used to construct well-known solutions and how it suggests the existence of 
novel, unstable solutions. We also discuss the possibihty of the existence of approxi- 
mate breathers (long-lived, spatially localized, temporally periodic configurations) in 
Minkowski spacetime. 



2.1 The Higgs-gauge Sector 

The bosonic sector of the electroweak theory is an SU{2) x U{1) gauged Higgs model. 
The Abclian coupling constant {g') is known to be much smaller than the non-Abelian 
coupling constant (g). For simplicity and clarity we choose g' = and ignore the 
dynamics of the U{1) hypercharge gauge fields. The three SU{2) gauge bosons are 
denoted by 1^^. These may be expressed as a matrix valued field, using the group 
generators r", 

W, = W^^. (2.1) 
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There are two complex scalar fields that form the Higgs doublet 



(2.2) 



where the subscripts denote the electric charge (had the U(l) fields been included). 
The Higgs may be written as a 2 x 2 matrix field 

^° 'A. (2.3) 

V -0; 00 J 

Under a gauge transformation, U{x) G SU{2), the fields transform as 

W^^U (w^ + '-d^^ t/t , $ ^ [/$ . (2.4) 
The Higgs-gauge sector is defined by the action 



A 



--(tr (<l>t$) -2i;2)' 



(2.5) 



where the field-strength tensor and the covariant derivative are defined as follows: 

W^, = d^W,-d,W,-ig[W„W,] , 

D^'i> = {d^-zgW,)<l>, (2.6) 

and A is the Higgs self- interaction coupling constant. 

The gauge symmetry is spontaneously broken because the Higgs doublet has a 
non-zero vacuum expectation value. This results in the following particle spectrum: 
a single scalar Higgs particle with mass m^^^ = 2vy/X and three degenerate gauge 
bosons with mass m^-' = gv/\^. The superscript '(0)' denotes that the masses are 
at tree level, and later in Chap. El we will discuss one-loop quantum corrections. 
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2.2 Solutions in Euclidean Spacetime 

In Euclidean spacetime, the Higgs-gauge action is positive definite and is given by 



(2.7) 



Since the metric is 6fj,^, there is no difi^erence between upper and lower indices. Also, 
for notational simplicity, we use the same Greek letters to denote spacetime indices 
in both Euclidean and Minkowski space. 

The classical equations of motion, obtained by extremizing the action with respect 
to the Higgs and gauge fields, are 



D F = ^ 



tr -Iv^ 



(2.8) 



The covariant derivatives for the Higgs and gauge fields are 



(2.9) 



2.2.1 Vacuum Solutions 

The totally-trivial configuration, $ = = 0, is obviously a solution at which 

the action has a global minimum of zero. All configurations gauge-equivalent to the 
totally-trivial configuration are also global minimum solutions. We refer to these as 
vacuum configurations. Any U{x) e SU{2) completely specifies such solutions as 
pure-gauge configurations 
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(2.10) 
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2.2.2 Non- vacuum Solutions and Topology 

In addition to the vacuum solutions with zero action, there are several solutions 
corresponding to local minima and saddle-points of the action. There are various 
strategies that enable us to find such solutions. One way is to embed known solutions 
of simpler theories into the electroweak theory. For example, the kink solution in 
1+1 dimensional, real 0^ theory becomes a domain wall solution in the electroweak 
theory. Also, Nielsen-Olesen vortices ^1] in the Abelian Higgs model become a family 
of string solutions in the electroweak theory ^Hl- A second strategy is to use topology, 
and this has lead to the discovery of the instanton fr\ and the sphaleron j^l . 

We now present a general prescription, which uses the topology associated with 
maps into the gauge group of the theory, to find known solutions in the theory and 
motivate the existence of novel ones. The basic idea, due to Manton, is to construct 
non-contractible loops in configuration space, and the "top" of the tightest loop should 
correspond to a solution. This is how the well-known electroweak sphaleron was con- 
structed. Klinkhamer has motivated the existence of other solutions by constructing 
different non-contractible loops [4 . We generalize this procedure and demonstrate 
how to construct all possible non-contractible loops in a gauge theory, and list the 
solutions indicated by these. 

2.2.3 Topological Prescription 

A finite Euclidean-action configuration is pure gauge at spacetime infinity: 

= vU , = ^Ud^U^ , (2.11) 

where [/ is a map from the boundary of spacetime to the gauge group, SU{2). We 
allow the configuration to have trivial dimensions, in which case it has finite action 
per unit volume of the trivial dimensions, and the domain of U is the appropriate 
subspace of the boundary of spacetime. For example, a static configuration has time 
as a trivial dimension and in order to have finite energy (action per unit time), it 
must be pure gauge at spatial infinity. 
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Now the third and fourth homotopy groups of SU (2) are non-trivial: 

Xl^{SU{2)) = Z,Iii{SU{2))^Z2. (2.12) 

So each map from 6*3 into SU (2) belongs to a homotopic class labeled by an integer 
winding number, and it cannot be continuously deformed into any map in a distinct 
class. Similarly, each map from 5*4 into SU{2) belongs either to the trivial class 
(which contains the trivial map) or the non-trivial class. Consider any topologically 
non-trivial map into SU{2) from ^'3 or 84^. Identify a subspace of the domain with 
the boundary of spacetime spanned by the non-trivial dimensions. Any remaining 
coordinates in the domain are interpolation parameters that define a sequence of pure- 
gauge, asymptotic configurations. These asymptotic configurations may be smoothly 
continued into the bulk to obtain a sequence of configurations. The sequence becomes 
a loop when we restrict all configurations on the boundary of the interpolation space 
to be the totally-trivial configuration ($ = vl, = 0). If the loop can be shrunk to a 
point, then the topologically non-trivial map into the gauge group can be continuously 
deformed to the identity, which is impossible. The top of the tightest non-contractible 
loop, obtained by minimizing the action (per unit volume of any trivial dimensions) 
for each point on the loop, should be an unstable solution. 

We now present detailed examples of the above construction and list the solutions 
indicated when all possible non-contractible loops are considered. 

2.2.4 Winding 1 Map from S3 to SU{2) 

Let /92, ct denote the angular coordinates that describe a three-dimensional sphere, 
S3, with < /3j < TT and < a < 27r. The 5*3 can be embedded as a unit sphere in a 
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four- dimensional Euclidean space and is described by unit vectors 



cos Pi 
sin Pi sin P2 cos a 
sin Pi sin P2 sin a 
y sin Pi cos y 



(2.13) 



Then, the canonical winding 1 map from 5*3 to SU (2) is given by the one-to-one map 



U^^\pi,p2,a) =nol + iniTi. 



(2.14) 



Now we will identify different subspaces of the domain S3 with the boundary of space- 
time, in accordance with the topological prescription, to find non-vacuum solutions 
to the classical equations of motion. It turns out that all the solutions found using 
this map are well-known, and our method simply encapsulates them into a unified 
framework. However, as we shall see later, other non-trivial maps lead to novel solu- 
tions. 

The Weak Instanton 

The entire 5*3 is identified with the boundary of spacetime by choosing the unit 
vectors n to be the unit position vectors x in Euclidean spacetime. This example is 
an exception because there are no parameters remaining in the domain to construct 
a loop. The asymptotic configuration is determined by 



U{x) = XqI + IXiTi 



(2.15) 



using eq. 12.111 Now we extend the configuration to the interior of spacetime, without 
introducing any singularities, using the ansatz 



fH{r)vU, 



= fwir)-Ud,U^ 



(2.16) 
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where r denotes the radius vector in Euchdean spacetime. The radial functions fn, fw 
go from at r = (for regularity at the origin) to 1 as r goes to oo (for finite- 
action). This configuration cannot be continuously deformed to the totally-trivial 
configuration because that would imply that U can be continuously deformed to 1. 
This indicates the existence of a topologically stable solution. Indeed, in the absence 
of the Higgs fields, the choice 



gives a stable solution to the equations of motion for every choice of the width w. This 
is the well-known weak instanton pp. It mediates non-perturbative fermion number 
violating processes via tunneling j2]. The action has a local minimum of Sir'^/g'^ at 
the weak instanton. 

When the Higgs fields are included, their contribution to the action can be made 
smaller by scaling to smaller distances, till the action reaches the no-Higgs value of 
Sn^/g'^ for which the configuration is singular. (This can be easily understood on di- 
mensional grounds.) This brings us to a crucial caveat of the topological prescription. 
The configuration space is a non-compact manifold, and so the non-contractible loops 
may run off to infinity without giving a solution. The method is useful to the extent 
that it suggests the existence of a solution, determines its stability (or lack thereof) 
and points to the region in configuration space where we can search for the possible 
solution. But it does not guarantee the existence of the solution. 

The Weak Sphaleron 

Now we consider time-independent configurations. We identify an 5*2 subspace of the 
5*3 domain of the winding 1 map (in eq. I2.14|) with the boundary of space: 



fw{r) 



r 



.2 _|_ yj2 



(2.17) 



(2.18) 
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where 6, (p are the polar and azimuthal angles respectively that span the spatial 
boundary. The remaining coordinate, ^, is the interpolation parameter. The multi- 
plication by e~*^'^^ fixes Uq = U-^ = 1, without changing the winding of the map. So 
defines a loop of asymptotic configurations, which can be smoothly continued into 
the bulk of space for each ^, to obtain a loop of configurations. For example, 

m = v{fHir)U^ + il-fHir))cos'^i), 
W,iO = fw{r)-U^d,Ul (2.19) 

where r is the spatial radius and the radial functions /h, fw go from at r = to 1 
as r — > oo. Suppose we could continuously deform the loop so that for every ^ the 
configuration is totally-trivial. Then f/*^^^(^, 6*, could be continuously deformed into 
the trivial map, which is impossible. So the loop is non-contractible. The classical 
energy along the loop starts at for ,^ = 0, then increases to some maximum value 
and finally goes down again to at ,^ = tt, where the configuration is totally-trivial. 
The loop can be made tighter, but cannot be shrunk to a point. The configuration at 
which the energy has a minimax (the top of the tightest loop) would be an unstable, 
static solution. We find that the minimax of our loop is at ^ = 7r/2 and corresponds 
to the solution to the following equations: 

= Vw{l-fw){l-Vw)-\9WfH{l-fw). (2.20) 

Near r = 0, oc r, fw oc r^. As r — > cxd, (1 — f^) oc l/(f r)exp(— 2-\/At>r) and 
(1 — /vk) oc exp [—gvr I This solution in fact solves the full equations of motion 
in eq. 12.81 and is the well-known weak sphaleron 01 C^- It is the lowest barrier 
between topologically inequivalent vacua (see Sec. IH.H.ll for this interpretation) and 
its energy determines the rate of fermion number violating processes at temperatures 
comparable to the electroweak phase transition scale. The radial functions fn, fw 
may be found numerically. For A = the energy of the weak sphaleron is about 
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8.607rv/g', which corresponds to about 7.6 TeV when we choose the experimental 
values V = 177 GeV, g — 0.63. (For non-zero A the energy is higher.) 

Wc should point out that had we ignored the Higgs fields, the minimax configu- 
ration could be driven to lower energies by expanding the size of the configuration. 
This is because the pure Yang- Mills theory has no scale (classically) and so the energy 
must be inversely proportional to the width of the configuration. The Higgs sector 
provides the vev scale, which prevents the non-contractible loop from running off to 
infinity. 

Weak Strings 

Finally we consider static configurations with one trivial dimension (say z) and iden- 
tify an S\ subspace of the domain with the boundary of the x — y plane. For 
example, 

u^uM = ft/('Ha,6,o)lVw(ei,6,¥'), (2.21) 



where ip is the azimuthal angle parameterizing the planar boundary, and the remain- 
ing domain coordinates, ^1,^2, are interpolation parameters. They parametrize a 
2-sphere (or 2-loop) of pure-gauge configurations at planar infinity. Note that the 
multiplication by C/^^^(^i, ^2, 0) ^ ensures that the boundary of the square spanned 

by .^1, .^2 is mapped to the identity, without afi^ccting the winding of the map, thereby 
making the 2-paramctcr sequence of maps into a 2-sphere of maps. The continuation 
into the bulk may be chosen to be 

^(6,6) = ^;[/i/(r)%,^, + {l-/^(r)}(l-sineisine2)l] , 

1^,(6,6) = fw{r)-U^.,i.d,Ul^^, (2.22) 

where r is the planar radius \/x^ + y'^. The radial functions /j^, fw go from to 1 as 
r goes from to infinity. Now we have a non-contractible 2-sphere of configurations, 
indicating the existence of an unstable solution at the top of the tightest sphere. We 
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find that the configuration at = ^2 = 7r/2 has the following non-zero components: 



0o = ^/H(r)e-^^ly=^ 



(2.23) 



If the radial functions are chosen to satisfy 




r ar 



(Pfw 



(2.24) 



then we have one of the well-known W-string solutions jlSl E|- This is in fact a 
Nielsen-Olesen vortex of the Abelian Higgs model jMj, embedded in a U{1) subgroup 
of the SU{2) theory. The above construction uses the U{1) subgroup generated by 
Ts, and other choices of the map U^i,^2iv) result in the solution embedded in other 
f/(l) subgroups, leading to a family of string solutions. As is the case for the weak 
sphaleron, the Higgs vacuum expectation value provides a scale that prevents the 
energy from approaching zero as the configuration width is increased to infinity. The 
instability of the weak string solutions detracts from their significance, especially 
with regard to electroweak baryogenesis where they could have played a crucial role. 
In Chap, m we discuss this in greater detail and explain how quantum effects could 
stabilize these configurations. 



2.2.5 Winding n Map from to SU{2) 

The above procedure can be repeated straightforwardly for a winding n{y^ 1) map 
from ^3 to SU{2y. 



If the entire domain is identified with the boundary of spacetime, then in analogy with 
the weak instanton, we should obtain a weak multi-instanton which is topologically 
stable and carries a topological charge of n. For time-independent configurations, 
we identify an S2 subspace of the domain with the boundary of space, with one 



f/(")(A,/?2,«)= \u^^\l3i,l32,a) 



1 



(2.25) 
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remaining interpolation parameter. This indicates the existence of a weak multi- 
sphaleron ^H] solution with one direction of instability. This would be the lowest 
energy barrier between two vacua that differ by winding number n. Finally, for 
static, planar configurations, we identify an 5*1 subspace of the domain with the 
planar boundary and obtain vorticity n W-strings with 2 directions of instability. 

We do not pursue these solutions further because they do not produce qualitatively 
different physical effects from those of the solutions obtained from the winding 1 
map. For example, the n-instanton would describe the production of n fermions 
via quantum tunneling through the sphaleron barrier. However, (n times) repeated 
instances of the 1-instanton already allows for this process, albeit with a possibly 
different probability amplitude. 

2.2.6 Non-trivial Map from ^4 to SU{2) 

Now we come to the relatively unexplored realm of classical solutions indicated by the 
non-trivial topology of maps from S'4 to SU{2). This topology is unique to SU{2), 
and doesn't exist for in the QCD gauge group SU{?>). Let (3i, P2, PsyCt denote the 
angular coordinates that describe a four- dimensional sphere, S'4, with < (3i < n and 
< a < 2tt. The map 

= e*'^-^[/W(/32,/33,a)e-^^-« [u^'\p„ (2.26) 

cannot be continuously deformed to the trivial map U = 1, and it belongs to the 
non-trivial homotopy class (as denoted by the superscript 'NT'). (Recall that U^^^ is 
a winding 1 map from S3 to SU{2), an example of which is in eq. 12.141 ) 

The /* 

Proceeding as before, when an 6*3 subspace of the S'4 domain is identified with the 
boundary of spacetime, we obtain a non-contractible loop of configurations in Eu- 
clidean spacetime. We ignore Higgs fields because on dimensional grounds, their 
contribution to the action approaches as the configuration is driven to a width 
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singularity. In the pure SU{2) theory, one example of the loop is 

W,{i) = f^{r)-U^d,Ul (2.27) 

where r denotes the radius in Euclidean spacetime. The action along the loop starts 
at at the totally-trivial configuration PF^(O), rises to a maximum and then returns 
to back at the totally-trivial configuration iy^(7r). The top of the tightest loop 
should correspond to a topologically trivial, unstable solution, the /* [1^]. It could 
contribute significantly to the path integral, since it extremizes the action. However, 
within the various ansatze that we consider, we find that the minimax corresponds 
to a widely separated instanton and anti-instanton pair. We have not succeeded in 
finding a new, localized solution. 

Another perspective on the possible /* solution comes from restricting fields to 
be totally-trivial at r — > cxd. Then the spacetime manifold TZ^ becomes compacti- 
fied to 5*4 for the purpose of maps into configurations. Vacuum configurations are 
pure-gauge and fall into two distinct homotopy classes, because Iii^{SU{2)) = Z2. 
Now, any interpolation from a vacuum configuration in the trivial class to a vacuum 
configuration in the non-trivial class must leave the vacuum manifold (otherwise the 
two configurations would be deformable into each other). The configuration with the 
smallest maximum action along all such possible paths would be a saddle point of the 
action. It would be the lowest action barrier between topologically inequivalent vac- 
uum configurations in Euclidean spacetime. This is analogous to the weak sphaleron, 
but in one higher dimension. 

The S* 

For static solutions we identify an 5*2 subspace of the 5*4 domain with the boundary 
of space and obtain a non-contractible 2-sphere of configurations. For example, 

u^„Uo,v) = [f/(^^)(6,e2,o,o)]V(^^)(ei,e2,^,¥^), 
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$(^1,6) = t;[/H(r)%,5, + {l-/^,(r)}(l-sineisin6)l] , 

= fwir)'-U^,,^,d,Ul^^, (2.28) 

where r, 6', ip denote spatial spherical coordinates. The radial functions fn, fw vanish 
at the origin and approach 1 as r ^ oo. The top of the tightest 2-sphere should give 
a static solution with 2 directions of instability, the S* |2Uj . 

However, we have not succeeded in finding a solution. As we tighten our 2- 
sphere of configurations within different ansatze, we find that the maximum action 
configuration approaches a separated sphaleron anti-sphaleron pair. 

The S*, if it exists, is a promising candidate to form a fermionic soliton that allows 
a heavy fermion doublet to decouple from the Standard Model (see Sec. 13. 5|) . 

The W — string* 

Finally, to obtain non-trivial solutions on the x — y plane, we can identify an Si sub- 
space of the 5*4 domain with the planar boundary. This would give a non-contractible 
3-sphere of configurations, which indicates the existence of a static, finite energy per 
unit length solution, with 3 directions of instability. Daunted by our inability to find 
the /* and the S*, we have not pursued this possibility. 

2.2.7 Caveats and Conclusions 

The topological prescription to find non-vacuum solutions to the classical equations 
of motion generalizes to any gauge theory. Once we know the gauge group, for each 
topologically non-trivial map into the gauge group, we construct a non-contractible 
n-sphere of configurations. One point on the sphere is chosen to be the totally-trivial 
vacuum configuration. Since the sphere cannot be shrunk to a point, the "tightest" 
sphere obtained by minimizing the action for each configuration on the sphere, could 
give solutions (corresponding to saddle points and local minima of the action). In the 
case of the simplified version of the electroweak theory, we have seen that this method 
allows us to find the well-known weak instanton, weak sphaleron, W-strings, and 
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their higher winding generahzations. It also suggests the existence of novel solutions, 
which are analogous to the well-known solutions, except constructed using a different 
topology. However, we have been unable to find these solutions. 

These topological arguments do not guarantee the existence of the solutions de- 
scribed above, because the configuration space is a non-compact manifold and the 
non-contractible loops may run off to infinity. For example, in the absence of the 
Higgs fields, the weak sphaleron's energy can be lowered by scaling to larger dis- 
tances and approaches as the solution approaches an infinite width singularity. 
Also, it not clear that two non-contractible loops obtained using different maps into 
the gauge group, give two distinct solutions. For example, we find that in our search 
for an /* using Il4{SU{2)) = Z2, the minimax configuration tends to break up into 
two instantons that were already discovered using Il3{SU{2)) = Z. Furthermore, the 
method is ignorant of solutions that have trivial gauge fields (such as the kink domain 
wall). Nevertheless, the topology points to possible solutions in the vast configura- 
tion space and once we know where to look, we can verify whether a solution exists. 
Furthermore, the method encodes whether the suggested solution is stable or not, 
and in the latter case it gives us the directions of instability. 

We find that all known and hinted time-independent solutions are unstable. The 
classical Higgs-gauge sector seems to have only sphalerons and no solitons. Of course, 
the existence of non-topological solitons (corresponding to local minima of the energy 
functional) cannot be excluded. However, having exhausted the topological properties 
of the theory, we are left with no guiding principle to enable a search for such objects. 
But if we consider quantum effects on the classical bosonic sector, then there are 
compelling reasons to expect the existence of quantum solitons and well-understood 
mechanisms to guide the search for them (see Chap. EI). 

2.3 Solutions in Minkowski Spacetime 

In this section we argue for the existence of approximate breathers in the electroweak 
theory. These are unnaturally long-lived, spatially localized configurations that are 
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periodic in real time. 



2.3.1 0^ in 1+1 Dimensions 

As a lead up to the electroweak breathers, consider the real 0^ theory in 1+1 di- 
mensions. This is arguably the simplest of all field theories, and yet it has enough 
structure to shed light on the mechanism that gives rise to breather configurations in 
classical field theories. The action is 

= / d'x ld^<l>d,<l> - ^ (0' - . (2.29) 

The equation of motion obtained by extremizing the action is 

- 0" = -A -v^)(f). (2.30) 

In addition to the vacuum solutions 

(j){x,t) = ±v, (2.31) 

this theory has the well-known solitonic kink (anti-kink) solutions 

(f)(x, t) = ±v tanh " , (2.32) 

^2(1 - 

centered initially at Xq and moving with velocity u. Now we shall see that there also 
exist approximate breathers in the theory. 

Consider small oscillations around the vacuum v, 

(l)^v + 5(f). (2.33) 

The linearized equation of motion is 

S4> - 6(1)" + 2Xv^6(l) = . (2.34) 
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Fourier transforming to momentum space, we get the hnear dispersion relation 



ujn„ = Vm^~+¥ , (2.35) 

where m = v^2Af is the mass of the scalar particles in the theory. The crucial 
observation is that there is a mass gap and the linear spectrum starts at m. 

Now we come to the second critical ingredient in the theory that allows for approx- 
imate breathers: the scalar potential is non-linear. More specifically, for increasing 
deviations of (p from v toward zero (the peak of the so called double- well potential), 
the curvature of the potential decreases. If this were a single particle potential, then 
this would imply that large amplitude oscillations around the vev will have a fre- 
quency lower than the linear frequency associated with the quadratic potential. This 
generalizes in the case of the field theory, and there are several spatially localized 
configurations which oscillate in time (or more precisely, the field value at the cen- 
ter oscillates in time) with a fundamental frequency below m. So, these oscillations 
should be stable against decay by linear modes radiation (i.e. boson radiation). How- 
ever, non-linearity also provides a flip side to this argument: the higher harmonics 
of the fundamental frequency must be present and they will fall in the band of the 
linear spectrum and cause the configuration to eventually decay. In discrete systems, 
the linear spectrum has an upper bound and the harmonics could be constructed to 
lie beyond the linear spectrum giving rise to stable discrete breathers. See |6j for 
an elementary review of discrete breathers, wherein these ideas are explained in the 
context of discrete systems. Also see [21] in which Gleiser has demonstrated the exis- 
tence of approximate breathers in continuum scalar field theories in 3-1-1 dimensions 
(for both symmetric and asymmetric double- well potentials). 

We now give an example of an approximate breather in the 1-1-1 dimensional 0^ 
theory. We set the dimensionless vev f to 1 . We also choose A = 1 and all dimensionful 
quantities are measured in units of VX. Consider the initial configuration 

0(x,O) = tanh2(x/w;) , (2.36) 
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with the width chosen to be w = 2. We start the above configuration at rest and then 
numerically evolve according to eq. I2.3(J1 Since our initial configuration and its time 
derivative are even functions of x, they remain even throughout the time evolution 
according to the equation of motion. So we consider only the positive half-line with 
vanishing spatial derivative boundary condition at the origin (in higher dimensions 
this condition is required for regular configurations). The numerical method is the 
following. We consider the x- interval between and Lq, where Lq is much larger 
than the initial width of the configuration. We discretize time and space with lattice 
spacings ht, respectively. Then the discretized equation of motion in eq. 12.301 
accurate to second order in time and space, is 

= + - 20i,j) + 2(j)ij - (pij^i - ht\(f)ij{(f)lj - f 2) , (2.37) 

where the first subscript labels the spatial lattice points and the second labels the 
temporal lattice points (so that (pij = (j){ihx, jht)). Note that the spatial lattice 
points within the interval are labeled hj i = 0, . . . ,N = Lq/H^. The above equation 
is explicitly in the form that allows 0(x) at the next time step to be determined if it 
is known at the current and previous time steps. For the spatial lattice point i = 0, 
the equation requires the value at the point i = —1 and we impose the condition that 
(f>-i,j = which corresponds to the function being even. At the wall (spatial lattice 
point A^), we impose the similar condition, ipN+ij = (/^n-ij- This condition ensures 
that no energy fiows past the wall (since it corresponds to vanishing derivative). 
To eliminate refiections from the wall returning to our local region of interest, we 
move the wall out at the speed of light (i.e. L = Lq + t) when we evolve an initial 
configuration. 

In Fig. l2-ll we plot the total energy, Ef, on the half-line as a function of time. This 
is conserved through the evolution. We also plot the energy, Ei, localized between 
a; = and x = 20. This falls very slowly and after a time of about 17,000, only 12% 
of the total energy has dissipated out of the local region. We do not yet understand 
what sets the scale for this decay rate. But it is clear that naturalness arguments are 



26 



violated in this example and the energy is localized for time periods much longer than 
all dimensionfull scales in the problem. In Fig. 12-21 we plot the value of the field at 
the origin as a function of time, towards the end of the time of evolution considered. 
This oscillates about cj) = v (asymmetrically because the potential is asymmetric 
about (f) = v). The fundamental period of oscillation is obtained by measuring the 
interval between peaks and the corresponding fundamental frequency turns out to be 
about 1.2. This is lower than the mass of the scalar particle, m = \/2 ^ 1.414, as 
expected. The existence of such approximate breathers seems to be a generic feature 
of the theory, and we have successfully constructed several such configurations with 
different frequencies and amplitudes. 
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Figure 2-1: The total energy, Et, on the positive half-line and the local energy, Ei, between 
X = and x = 20, as a function of time for the initial configuration in eq. 12.361 in the (j)'^ 
theory. All dimensionful quantities are in units of \/A and the vev u is 1. 

We extend the 0^ theory by considering complex and gauging the global U{1) 
symmetry. The gauge boson becomes massive after eating the Goldstone degree of 
freedom. So the linear spectrum has mass gaps corresponding to the masses of the 
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Figure 2-2: The field value at the center of the configuration as a function of time, for the 
initial configuration in eq. 12.361 in the cj)'^ theory. All dimensionful quantities are in units of 
\/A and the vev u is 1. 

scalar and vector bosons. We find that the breathers in the real (p^ theory are stable 
against perturbations in the complex gauged theory, as long as the U{1) charge is 
large enough so that the fundamental breather frequencies are lower than the mass of 
the gauge boson (which is proportional to the charge). However, when we reduce the 
charge so that the fundamental breather frequency falls within the spectrum of small 
gauge field oscillations, the breathers dissipate their energies rapidly. This is further 
evidence that the existence of approximate breathers is due to the presence of a mass 
gap and a non-linear potential that allows oscillations with a frequency smaller than 
the mass of the lightest particle in the theory. 
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2.3.2 Breathers in the Electroweak Theory 

Now that we understand the mechanism that allows long-lived, spatially localized, 
temporally periodic configurations to exist in a classical field theory, we can ask if 
there is any possibility for such objects to exist in the Standard Model. The answer is 
yes. When we ignore the hypercharge gauge fields in the electroweak theory and con- 
sider the SU (2) Higgs theory (as described in Sec. l2.1|) . there is no unbroken subgroup, 
and all three W gauge bosons are massive. So, the theory fulfills the first require- 
ment of having a mass gap in the linear spectrum. Also, the Higgs potential is the 
usual sombrero potential and so it seems likely that we could set up large-amplitude 
configurations with oscillation frequencies in the mass gap. So the theory appears to 
have all the ingredients required for the existence of approximate breathers. When 
we include the U{1) sector, the photons remain massless after symmetry breaking. 
However, it is conceivable that the breathers can be made electrically neutral so that 
they don't dissipate by electromagnetic radiation. 

We are currently investigating whether such objects exist in the electroweak the- 
ory. They could have many significant implications. Firstly, since the lifetimes of the 
breathers is expected to be many orders of magnitude larger than all natural scales in 
the problem, they could shed light on naturalness in field theories. Secondly, if these 
breathers could exist just after the electroweak phase transition, they would set up 
out-of-equilibrium regions in space, with implications for electroweak baryogenesis. 
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Chapter 3 
Quantum Solitons 



We explain how certain Higgs-gauge configurations may be stabilized in the Standard 
Model, by carrying heavy quark quantum numbers. These fermionic solitons would 
maintain anomaly cancellation in the low-energy effective theory and allow heavy 
fermions to decouple from the chiral gauge theory. We describe the methodology of 
looking for such objects in the configuration space. We present a technique based on 
scattering theory that enables an efficient, exact and properly renormalized calcula- 
tion of the one-loop effective energy, which is one crucial ingredient in the stability 
analysis of quantum solitons. Finally, we describe our search for a fermionic soliton 
within a spherical ansatz. Much of this investigation originally appeared in jU]. 

3.1 The Idea 

A topological soliton is a non-vacuum, static configuration that is topologically sta- 
ble. It carries a conserved topological charge which prevents it from decaying into a 
vacuum configuration with no topological charge. Analogously, it may be possible to 
stabilize a configuration by making it carry a conserved quantum number. We use the 
term quantum soliton to refer to any such quantum-stabilized object. In this section 
we briefly describe the mechanism and motivation for these solitons, and elaborate 
on these ideas in the following sections. 

There is a natural mechanism in the electroweak theory for the existence of quan- 
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turn solitons stabilized by carrying fermion number (say top quark number), which 
stems from the polarization of the fermionic vacuum in a background of gauge and 
Higgs fields. Several configurations are known to tightly bind fermions in their vicin- 
ity. These configurations consist of classical solutions (the sphalerons discussed in 
Sec. I2.2.2|) as well as non-solutions. The existence of tightly-bound levels suggests 
that it may be energetically favorable for a certain number of fermions (say Nf) to be 
trapped by such backgrounds, with a small associated occupation energy, Ei^/\ The 
binding energy could outweigh the cost in classical energy, Ed, to set up the config- 
uration, i.e. Eel + Eocc'' < rufNf, where m/ is the mass of the fermion. However, to 
be consistent to order h, we must also include the fermion Casimir energy: the renor- 
malized energy shifts of all the other fermion modes. For static configurations this 
is the renormalized one-loop fermion vacuum energy -Evac- Thus, the minimum total 
energy of an arbitrary configuration C = {$, Wf^}, which carries fermion number Nf, 
is 

eS'\c] = E4C] + Ei^f^lC] + i?.ac[q . (3.1) 

The subscript 'eff' denotes that the energy is an effective energy obtained by in- 
tegrating out the fermions from the theory, as explained in Sec. 13.21 If C decays 
into a vacuum configuration then it must create Nf perturbative fermions (ignoring 
anomalous violation of fermion number, which is exponentially suppressed). How- 
ever, the strong fermion binding suggests that there exist configurations such that 
El^f\c] < rrifNf. Then, C can decay only into a quantum soliton at which E^^f^ 
has a local minimum. 

The existence of such fermionic quantum solitons would provide an attractive 
resolution to the decoupling puzzle in the Standard Model (and other chiral gauge 
theories) [22] ■ A fermion obtains its mass through Yukawa coupling to the scalar 
Higgs via the well-known Higgs mechanism. Explicit mass terms are prohibited by 
gauge-invariance. So, as we increase the mass of a fermion (thereby making the 
denominator in the propagator suppress loop corrections) we also increase the Yukawa 
coupling which gives a corresponding enhancement from the vertices. Moreover, the 
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heavy fermion cannot simply disappear from the spectrum because then anomaly 
cancellation would be ruined. However, it is plausible that the large Yukawa coupling 
gives rise to a quantum soliton in the low energy theory. This carries the quantum 
numbers of the decoupled fermion and maintains anomaly cancellation using the 
mechanism described by D'Hoker and Farhi in the SU{2)l Skyrme model [23] ■ 



3.2 The Higgs-gauge Sector Effective Energy 

In addition to the bosonic sector described in Sec. 12.11 the Standard Model has a 
rich fermion sector of three generations of quarks and leptons. We will be primarily 
interested in one-loop quantum effects of heavy fermions on the Higgs-gauge sector. So 
we ignore inter-generation mixing and set the Cabibbo-Kobayashi-Maskawa (CKM) 
matrix to the identity. Let \1/ denote an isospin doublet, say the top-bottom quark 
doublet. Only the left-handed fermions are charged under a gauge transformation, 

u{x) G su{2y. 

^L^U^L,'^R-^'^R, (3.2) 

where the subscript L, R denote left and right handed projections respectively. The 
action for the doublet is 

Sf[^, $, W^] = J d^x [^Lt^^D^^L + ^Ri^d^^R - f {^L^^R + h.c.)] . (3.3) 

The covariant derivative is 



D^^L = {d^-tgW,)^L, (3.4) 

and / is the Yukawa coupling constant (which may be different for each doublet). 
For simplicity we have used the same Yukawa coupling for both components of the 
doublet, which results in the degenerate mass rrif = fv, when the Higgs develops 
the vev v after spontaneous symmetry breaking. In the absence of inter-generation 
mixing, the fermion sector Lagrangian is given by the sum of the Lagrangians for 
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each fermion doublet. We introduce the potential 

V{<l>, W^) = -gj>^W^{x)^^^ + f{h{x) + tvp''{x)T-^5) , (3.5) 



where h{x),p°'{x) are real functions that characterize deviations of the scalar fields 
from the vev: 

<!>{x) -V = h{x) +ivp''{x)T'' . (3.6) 
We write the Lagrangian density as the sum of a free part and an interacting part: 

Sf[^, V{^, W^)] = J d^x {i-i^'d^ - m/) ^ - W^] . (3.7) 

Since we have ignored the U(l) hypercharge sector, there are no perturbative 
anomalies associated with triangle diagrams, and we only have Witten's global anomaly 
[21] to contend with. This requires an even number of SU {2)l fermions in the theory, 
otherwise the theory is not defined. (In the path integral formulation, the partition 
function becomes zero. Alternatively, in a Hamiltonian formulation, the theory has 
no gauge-invariant states |2S1-) For now we restrict our attention to a single, heavy 
fermion doublet and ignore the lighter spectator doublets present for anomaly can- 
cellation. We investigate if the heavy doublet can be decoupled from the theory by 
becoming solitonic at a lower energy scale. The full theory we consider is then given 
by the Higgs-gauge sector, Sh, described in Sec. 12.11 and the heavy fermion sector, 

W^, % = W,] + Sf[^, V{^, W,)] . (3.8) 

We integrate out the fermion doublet from the theory to obtain the effective action 
for the Higgs-gauge sector: 

The normalization has been chosen so that the effective action is equal to the classical 
action for vanishing interaction potential, V ^ defined in eq. fj3.7p . If iS^^ denotes the 
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Feynman diagram with one fermion loop and n external insertions of {—iV{(^, ^a*)); 
then 



n=l 



(3.10) 



The Feynman diagrams can be computed in a prescribed regularization scheme. The 
divergences that emerge as the regulator is removed are canceled by counterterms. 
We introduce the renormalized and counterterm actions, S'f^^^ and S'j^*'* respectively, 
by expressing the bare parameters in Sh in terms of renormalized parameters and 
counterterm coefficients to obtain Sh = S'f^'^^ + S'ff^. The renormalized action is 
the original Higgs-gauge sector action, eq. ()2.5|) . with renormalized parameters sub- 
stituted. It can be written as the sum of all possible renormalizable, gauge-invariant 
terms: 



S 



(ct) 
H 



+C3 



citr iW'W^,) + C2tr 

2' 



+ C4tr 



(3.11) 



The coefficients q depend on the regulator. For notational simplicity we do not 
introduce a different notation for the renormalized parameters and fields ($, W^^ g, A). 

If we consider static Higgs-gauge fields (in the Wq = gauge) and restrict time to 
the interval T, then the effective energy functional is 



1 — *^oo ± 



(3.12) 



where Ec\ refers to the classical energy of the Higgs-gauge sector: 



Ai[$, Wi\ = \ d'xl -tr {W,,W,,) + -tr ([A$]' A<f j + - 



The fermionic vacuum energy is 



tr $T$ _ 2v 



(3.13) 



Aac[$, Wi\ = 4"d['^', + ^ct[$, W,] 

n=l 



(3.14) 
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where each regulated Feynman diagram contribution is 



(3.15) 



and the counterterm contribution is 



Eet[$, W,] - j (i^xj-citr (Wi^Wi^) + C2tr ([A*]^ A*) 



-C3 



tr ($1"$) - 2v^] - C4tr [$1"$ 



(3.16) 



The entire one-loop effective energy receives contributions also from gauge and Higgs 
loops. We ignore these contributions because we believe the fermion loops are fun- 
damentally responsible for the phenomena associated with fermion decoupling. If we 
imagine that the fermions have Nq internal degrees of freedom (e.g. color), then this 
approximation becomes exact for large Nc- Nevertheless we set Nc = 1. 



3.2.1 The Counterterms 

The counterterms render finite by canceling the divergences in E^, for n = 1 
through n — 4, that arise when the regulator is removed. To unambiguously de- 
termine the counterterm coefficients, q, we impose physical on-shell renormalization 
conditions: 

a. We choose the vacuum expectation value (vev) of h{x) to be 0, which ensures 
that the vev ($) = vl stays fixed at its classical value and the perturbative 
fermion mass does not get renormalized. This is equivalent to a "no-tadpole" 
renormalization condition and determines C3. 

b. We fix the pole of the Higgs propagator to be at the tree level mass, mh = rii'h^ , 
with residue 1. These conditions determine C2 and C4. 

c. We have various choices to fix the remaining undetermined counterterm coeffi- 
cient Ci. We choose to set the residue of the pole of the gauge field propagator 
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to 1 in unitary gauge. Then the position of that pole, i.e. the mass of the gauge 
field, prediction. 

The resulting counterterm coefficients, listed in Appendix El As explained 

under item c, the mass of gauge fields is constrained by the other model parameters 
when fermion loops are included. With our choice of renormalization conditions, it is 
the solution to the implicit equation 



ml 



1 + rdxxHi-x?^^ 



+6 / dxx{\-x)\Y^ ^ \ \ dxln 

JO mj- JO 



A(x, 
m) 



, (3.17) 



with A(x, g^) = — x(l — x)g^. Recall that m^-* = gvj \/2 is the tree level pertur- 
bative mass of the gauge fields. 

3.2.2 The Fermion Vacuum Energy 

We briefiy summarize methods introduced in refs. f7| IB] (see ref. ^] for a review 
and a list of additional references) that enable an exact calculation of one-loop ef- 
fective energies in quantum field theories. We use these techniques to compute the 
fermion quantum corrections (-Evac) to the energies of Higgs-gauge configurations. As 
an aside, we want to point out that we have fruitfully used these methods in the study 
of vacuum energies of quantum fields in the presence of boundary conditions, and the 
resulting Casimir forces and stresses on the boundaries j^Hl HZl I2H] • The boundary 
conditions give rise to additional divergences in such calculations, and renormaliza- 
tion is a contentious issue. We replace the boundaries with background fields that 
couple to the quantum fields in such a way that in an appropriate limit, they con- 
strain the quantum fields to obey the boundary conditions. Then we calculate in 
this renormalizable field theory to obtain properly renormalized, finite results in the 
presence of the background. Finally, we take the boundary condition limit on the 
background to obtain the Casimir energy in the presence of boundaries. We find that 
the Casimir energy diverges in the boundary condition limit and this divergence can- 
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not be absorbed into a renormalization of the parameters of the theory. This means 
that the Casimir energy and other quantities such as the stress on an isolated surface 
are sensitive to the physical ultra-violet cutoff. On the other hand, the Casimir force 
between objects and the energy density away from the surfaces are finite and well- 
defined. Since these topics lie outside the realm of the electroweak theory, we have 
omitted a more detailed discussion of this here. 

In Eq. fj3.14|) . the fermion vacuum energy in a Higgs-gauge background is for- 
mally expressed as an infinite sum of Feynman diagrams. We are interested in heavy 
fermions with large Yukawa couplings, which forbids a perturbative calculation in the 
couplings. Also, we will usually consider configurations with widths of the order of 
the fermion Compton wavelength, which does not allow an approximate calculation 
using a derivative expansion. In order to compute non-perturbatively, we make use 
of the fact that E'vac is also given by a sum over the shift in the zero-point energies 
of the fermion modes due to the background fields. We write this formal quantity as 
a sum over bound state energies, ej, (times their degeneracies, Dj) and a momentum 
integral of the energy of the continuum states weighted by the change in the density 
of continuum states, Ap(/c), that is induced by the background fields, 

1 1 f°° I 

E..c = -^J2Dj\^j\-2i dk^k^ + m}Apik) + E,,. (3.18) 

The momentum integral and E^t are both divergent, but their sum is finite because 
the theory is renormalizable. We render the integral finite by subtracting the first 
terms in the Born series expansion of the density of states and adding back in exactly 
the same quantity in the form of Feynman diagrams: 

^vac = -^E^jl^^-^i dk^k^ + m}^Apik)-J2Ap^Hk) 

N 

+ J2E^l + E,,. (3.19) 

1=1 

When we combine E^ with J2iLi -E'fd) cancel the divergences as well as implement 



37 



the renormalization prescription. As a result, the above expression is manifestly finite. 
The minimal number of required Born subtractions, N, is easily determined by an 
analysis of the superficial degree of divergence of the one-loop Feynman diagrams. 
For our theory, = 4. 

We will work with background fields in the spherical ansatz [221, as described 
in Sec. 13.4.11 Then we can express Ap{k) (and its Born series) in terms of the 
momentum derivative of the phase shifts j^, induced by the background fields, of 
the Dirac wave-functions. 

Here we have expanded the S'-matrix in grand spin channels labeled by G. The 
grand spin is the vector sum of total angular momentum and isospin. It is conserved 
by the potential, eq. ()3.7j) . for spherically symmetric field configurations. Note that 
this method requires us to consider sufficiently symmetric configurations so that the 
scattering matrix can be expanded in partial waves, and we have chosen spherical 
symmetry. We will later consider cylindrical symmetry when we investigate elec- 
troweak strings in Chap. 01 The asymptotic scattering states are labeled by parity 
(— l)*^ and total spin G ± 1/2 and we obtain a 4 x 4 S'-matrix in general (except in 
the G = channel, where it is 2 x 2). We let Sc^aik) denote the sum of the eigen- 
phase shifts at momentum k in channel G and a = ± specify the sign of the energy 
eigenvalue, u = ±^/c^ + m'j. Note that the single particle spectrum is not symmetric 
because the Dirac Hamiltonian is not charge conjugation invariant. The degeneracy 
is given by Dq = 2G + 1. In Appendix^ we show in detail how to use the Dirac 
equation to calculate the bound state energies and scattering phase shifts needed in 
the computation of -Evac- 

To simplify the calculation we only subtract the first two Born approximants and 
eliminate the remaining log-divergence in the momentum integral by using a limiting 
function approach developed in ref . [SI] . The final expression for the fermion vacuum 
energy is (after integrating by parts and using Levinson's theorem which relates the 
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phase shifts at threshold to the number of bound states) 

^vac = -^E(2G, + l)(|e,|-m;) + i?(i'2)+E(3'^) 
i 

+ r ^ . ~m + ^SUO) , (3.21) 

where Gj is the grand spin associated with the bound state j and 

OD 

m= E E(2G + l)(5c,.(A:)-C(^)-CW)+W^)- (3.22) 

(7e{+ ,-} G=o 

Here 6^Q^{k) denotes the z**^-term in the Born series of Scaik)- After subtracting 5^^^ 
and 5^"^^ from 5, the momentum integral does not contain any contributions that are 
linear or quadratic in V. The limiting function for the sum over all eigenphase shifts, 
'^iim(^), eliminates the logarithmically divergent pieces that are third and fourth order 
in V from the momentum integral. Its analytic expression can be extracted from the 
divergent pieces of the corresponding Feynman diagrams and is given in Appendix 
1X1 ■ eq. flA.22j) . Furthermore, E'^^''^^ denotes the contribution up to second order in 
V from the renormalized Feynman diagrams. Its explicit expression is displayed in 
Appendix^ Finally, E^^'^^ contains the third and fourth order counterterm contribu- 
tion combined with the divergences in the third and fourth order Feynman diagrams 
that have been subtracted from the momentum integral via (5iim. Again, its explicit 
form can be found in Appendix ^ 

Thus we compute an exact, finite, renormalized, gauge-invariant effective energy 
functional, £'eff[$,Wj], defined in eq. (|3.12p . for the Higgs-gauge sector, with the 
fermion fields integrated out. 



3.3 The Energy of a Fermionic Configuration 

We are interested in exploring the possibility of the emergence of a stable, fermionic 
soliton in the Higgs-gauge sector of the theory as we increase the Yukawa coupling 
(thereby making the perturbative fermion heavier). In the previous section we out- 
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lined the procedure that allows us to calculate the effective energy when the fermions 
are integrated out. Now we analyze the minimum additional energy required to as- 
sociate unit fermion number with a particular Higgs-gauge configuration C, where 
C = {^,Wi}. First in section 15.3.11 we determine the integer fermion number F[C] 
of a configuration. (This is subtle because we have to contend with the anomalous 
violation of fermion number). Then we occupy or empty levels of the single-particle 
Dirac Hamiltonian to obtain the lowest energy state with net fermion number 1. If 
F[C] = 0, then the lightest positive bound state needs to be filled and the occupa- 
tion energy E^]l = eiowest, where the superscript '(!)' denotes that levels have been 
occupied/emptied to obtain fermion number 1. If F[C] = 1 then Ej^H = because C 
is already fermionic, and so on. Thus, the minimum total energy of a single fermion 
associated with a configuration is 

= E,,[C] + E^^^iC] + EilliC] . (3.23) 

In previous works, such as the fermion vacuum contribution was omitted 
from the above equation and stable non-topological solitons were found. We consider 
such calculations inconsistent, because E^H and -Evac are both order h. We will see 
explicitly in section 13.4.31 that E^^c makes a significant positive contribution when the 
Yukawa coupling is large enough that the perturbative fermion mass is comparable 
to the classical energy. 

Since we refer to all these different energies frequently in the rest of the paper, we 
summarize our notation in Table IXTl 



Ed 


Classical Higgs and gauge energv. eq. ()3.13|1 


E 


Renormalized fermion vacuum energy, eqs. (13.141 13.181) 


EeS 


One-fermion-loop effective energy, Ed + -Evac 


occ 


Smallest occupation energv for fermion number m. eq. ()3.23|) 




Smallest effective energy in the fermion number m sector, E^s + -E'occ'* 



Table 3.1: Definitions of some of the energies which appear in our analysis. 
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3.3.1 The Fermion Number of a Configuration 

First we review properties of the Higgs-gauge configuration space and the classical 
energy functional defined on it. From the expression for Ed in eq. ()3.13|) . it follows 
that configurations 

Wj = V(")9,f/(")^ (3.24) 

have Ed = and we refer to them as vacuum configurations. Here U^^^ is any map 
from to SU (2) with winding number n. We have compactified position space from 
7^.3 to 5*3 by restricting the fields to be totally-trivial at spatial infinity. Vacuum 
configurations with the same winding number are equivalent under small (winding 
number 0) gauge transformations. We use C-"^ to denote the homotopic class of 
vacuum configurations with winding number n. Topologically inequivalent vacuum 
configurations are related by large (nonzero winding number) gauge transformations. 
Along any continuous interpolation between two configurations, one in C(") and the 
other in C*^™-* (with n 7^ m), there exists a configuration, C, with maximum classi- 
cal energy. Since no f/*-"-* can be continuously deformed into any U^"^\ Ed[C] > 0. 
The configuration corresponding to the minimax of these energies, when all interpola- 
tions are considered, is the classical sphaleron. This is the weak sphaleron discussed 
in Sec. 12.2.41 where it was constructed using the topological prescription of non- 
contractible loops. When the fermion vacuum energy is added to the classical energy 
to obtain the effective energy {Ees = Ed + -Evac), not only does the magnitude of 
the minimax energy change, but its location in configuration space shifts as well. We 
therefore define the quantum- corrected sphaleron to be the configuration that has the 
lowest of the maximum effective energies along all interpolations between topologi- 
cally inequivalent vacuum configurations. 

We associate any configuration C with a unique class of vacuum configurations 
by continuously deforming C in the direction of the negative gradient of the classical 
energy until we get a configuration in C^'^^ for some n = n(C). We call C^"*^*"-*^ the 
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connected C-class of C and we say that C is in the winding number n basin. Note that 
the classical sphaleron and all configurations that descend to it are on the boundary 
between different basins. Also note that we do not fix the gauge during the gradient 
descent, otherwise all configurations would descend to the same topological vacuum. 

For any two configurations Ci and C2, we define the spectral flow S[Ci,C2\ to be 
the number of eigenvalues (levels) of the single particle Dirac Hamiltonian that cross 
zero from above minus the number that cross zero from below along any interpola- 
tion from Ci to C2. See jSS] for an early discussion on the relation between level 
crossings and the anomalous violation of fermion number. The fermion number of a 
configuration C is defined as 

F[C] = S[Ci,C] (3.25) 

with Ci E C^'^^'"^\ Since the Dirac spectrum is gauge- invariant, F[C] does not de- 
pend on which particular Ci is chosen from the connected C-class. Moreover, F[C] is 
gauge-invariant even under large gauge transformations. Also, F[C] does not depend 
on the chosen interpolation because the spectral fiow is the same for all interpola- 
tions. This definition of the fermion number can be readily understood for a C that 
has classical energy less than the energy of the classical sphaleron. A continuous 
interpolation from any configuration in C'-"*^'"^-' to C preserves net fermion number 
because anomalous fermion number violations require the Higgs-gauge fields to cross 
the sphaleron barrier. However, this does not preclude level crossings corresponding 
to pair production during the interpolation. Thus, defining vacuum configurations 
to have zero fermion number leads to eq. ()3.25|) . Configurations that have classical 
energies larger than the classical sphaleron are not separated by an energy barrier 
from topologically inequivalent basins, so it is not clear what their fermion number 
should be, although our definition does assign a unique F to them. 

Having determined the fermion number of a configuration, we can use eq. ()3.23p 
to find E^^, the minimum effective energy of a configuration in the fermion number 
1 sector. 
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3.3.2 Stability of the Soliton 

We would like to know if there exists a configuration at which the one-fermion effective 
energy functional, E^^ , has a local minimum. This configuration would be a fermionic 
soliton. We carry out a variational search, looking for a configuration C such that 
-^eff [^] < ^^"i -^eff [^] < -^q.s.) whcrc i^q.g, is the effective energy of the quantum- 
corrected sphaleron. The first condition ensures that C cannot simply decay into 
a configuration with zero classical energy plus a perturbative fermion. The second 
condition prevents C from rolling over the quantum-corrected sphaleron, giving up 
its fermion number and then rolling down the E^g surface to a vacuum configuration. 
Finding a configuration with these properties would guarantee the existence of a 
nontrivial local minimum of E^V . 



3.4 The Search for a Spherical Soliton 

In this section we describe our search for the soliton. We first review the spherical 
ansatz for the gauge and Higgs fields. We then outline the restrictions imposed on the 
variational ansatze used to search for a soliton. Finally, we report on our search within 
two physically motivated sets of ansatze: the "twisted Higgs" and "paths over the 
sphaleron". Note that throughout this section the perturbative fermion mass is set 
to 1 so that energies and lengths are measured in units of and l/mf, respectively. 

3.4.1 The Spherical Ansatz 

We only consider static gauge and Higgs fields in the spherical ansatz. This enables us 
to expand the fermion S matrix in terms of partial waves labeled by the grand spin G, 
as described in Appendix ^ Our method for calculating the fermion vacuum energy 
requires such an expansion. Under these restrictions (and in the Wq = gauge, which 
for smooth fields is obtained by a non-singular gauge transformation), the fields can 
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be expressed in terms of five real functions of r: 



WJx) = -WHx) 



29 



. X . . a{r) , . . . 7(r) 



= V [s(r) + ip{r)TjXj] 



(3.26) 



where x is the unit three-vector in tlie radial direction. 

The ansatz transforms under a f/(l) subgroup of the full SU{2) gauge symmetry, 
corresponding to elements of the form 



g{x) = e*-^M-^-^^/2 ^ 



(3.27) 



with Oi transforming as a 1 dimensional vector field, s + ip as a. scalar with charge 
1/2, and a + ^(7 — 1) as a scalar with charge 1. It is convenient to introduce the 
moduli p, S and phases 9, r] for the charged scalars: 



- ipe''^ = a + i(7 - 1) and Se*'' = s + ip . 



(3.28) 



From now on we will specify a configuration using the five functions ai{r) , p{r) , 6{r), S(r) 
and ?7(r). Regularity of Wi{x) and $(x) at x = requires that 



P(0) = 1, 

P(0) = 0, 

^(0) = -2ng7T, 

am = 9'{0), 

either S(0) = or r]{0) = -n^vr . (3.29) 

Here no^Uj^ are integers and primes denote derivatives with respect to the radial 
coordinate. 

For the gauge transformation g{x) in eq. ()3.27|) to be non-singular, we require 
/(O) = —2n7c, where n is an integer and we denote this boundary condition as a 
superscript: /(r) = f^^\r). If f^'^\r) is restricted to be as r — 00 (which is 
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equivalent to g{r- — > oo) = 1) then n is the winding of the map g{x) : S'^ SU{2). 
So the topology of the vacuum configurations persists in the spherical ansatz. The 
classical energy of eq. ()3.13|) takes the form 



Eel 



2 tot 



aiY + 



2r2 



+ 



((p-l)^ + Vsin^^ 



and winding number n vacuum configurations of eq. ()3.24j) now become 



(3.30) 



p{r) = 1, S(r) = l, 

0{r) = /")(r), r^(r) = ^^, ai(r) = /(")'(r) . 



(3.31) 



We want the Higgs-gauge fields to have finite classical energy. So we require a field 
configuration of the form eq. ()3.3H) as r — oo, and the restriction that /'•"'^(oo) = 
uniquely specifies the boundary conditions on the fields at infinity. At x = 0, the 
boundary conditions on p (specified in eq. fl3.29|) ) make the energy density finite and 
we do not require any additional constraints. 

The anomalous violation of fermion number is given by the anomaly equation 



(3.32) 



where i^'^ is the Chern-Simons current. It is useful to consider the charge associated 
with it: 



9' 



cs 



2tt Jo 



I d\ii (WidjWk - '^tgWiWjWk 



ai + p'^{9' - ai] 



(3.33) 



This is a non-integer in general, and is equal to the integer winding number of Z*^"-* 
for the configurations of eq. 1)3.311) . and a half-integer for the sphaleron ^H]- Under 
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a winding n gauge transformation, A^cs ~^ ^cs + For background fields that 
interpolate between topologically distinct vacuum configurations, the net fermion 
number produced is given by the change in A^^cs. 

3.4.2 Restrictions on the Variational Ansatze 

Our methods allow us to consider any static, spherically symmetric configuration, C, 
in the Higgs-gauge sector, specified by the five real functions cii(r), p(r), 9{r), S(r) 
and ri{r). In principle, we could numerically minimize the fermionic energy, i?gg [C], 
in terms of the five functions and determine if a soliton exists. In practice however, 
that is not feasible. So instead we limit ourselves to the variation of a few parameters 
in ansatze motivated by physical considerations. 

We restrict our variational ansatze to those that obey the above described bound- 
ary conditions at the origin and at infinity. In addition, we restrict the Higgs fields 
to lie within the chiral circle, S(r) < 1, because otherwise the effective potential is 
unbounded from below. (The leading terms in the derivative expansion of eq. 13.101 
can be found in jHl])- Finally, the effective theory (obtained by integrating out the 
fermions) has a Landau pole in the ultraviolet, reflecting new dynamics at some cut- 
off energy scale or equivalently at a minimum distance scale. Conflgurations that are 
large compared to this distance scale are relatively insensitive to the new dynamics at 
the cutoff, but smaller conflgurations are sensitive. For small widths and large cou- 
plings, the Landau pole becomes signiflcant and leads to unphysical negative effective 
energies in eq. (I3.12p jSHl ES] ■ We have to be wary of this in estimating the reliability 
of our results. See for more detailed discussions on the effective potential and 
the Landau pole. 

3.4.3 Twisted Higgs 

We flrst consider twisted Higgs conflgurations, with = 1 so that rj^r) goes from — tt 
at r = to as r ^ oo. The other functions are trivial: ai(r) = 0, p(r) = 1, 6{r) = 
and S(r) = 1. If we smoothly interpolate from a vacuum conflguration (in the 
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connected C-class) to such a configuration, we observe that one fermion bound state, 
that originates in the positive continuum, has its energy decrease sharply. The wider 
the final twisted Higgs configuration, the closer this level ends up to the negative 
continuum. At a width of order l/rrif, it has energy zero, eliminating any occupation 
energy contribution, E^H, to the fermionic energy, E^^ , associated with it. The 
existence of this level makes such twisted Higgs configurations attractive candidates 
for the variational search. 

We consider one such twisted Higgs configuration with a width characterized by 
a variational parameter w, 

rj = -T.e-'^l'" (3.34) 

and add various perturbations to it. For instance, one among many of our variational 
ansatze (in the 6' = gauge) is a four parameter ansatz with parameters po? • • • iV^'- 

1 + [r/wY 
1 + [r/w) 
1 + [r jwy 

P = 1+P3 .1'('^f,3 e-/-, (3.35) 
1 + [r/wy 

where — 1 < pi < (to keep the Higgs field within the chiral circle and its magnitude 
positive) and p3 > —5.23 (to keep p positive). For a prescribed set of theory parame- 
ters {rriwimh and /) we determine the gauge coupling g = \/2'm^^'^ /v from the renor- 
malization constraint eq. (|3.17|) . We then vary the ansatz parameters (w,po) • • • ?P3) 
to lower the fermionic energy e'§. We find that the gain in binding energy is insuf- 
ficient to compensate for the increase in the effective energy E^.^. Through all our 
variations, E^^^ is strictly greater than m/ and we find no evidence for a soliton. The 
same result was obtained in [HI] without gauge fields, and the extra gauge degrees of 
freedom do not seem to help in the twisted Higgs ansatz. 

We also find that the fermion vacuum contribution, -Evac, to the fermionic energy. 
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E^g destabilizes would-be solitons. Consider a linear interpolation from the trivial 
vacuum configuration to the twisted Higgs configuration in eq. fl3.34|) with gauge fields 
set to zero. We introduce the interpolating parameter ^ which goes from to 1: 



We choose the Yukawa coupling to be / = 10 and the Higgs mass to be v/\/2. Since 
the gauge fields are trivial, the classical energy as well as the Dirac spectrum are 
independent of the gauge coupling g and the gauge bosons mass m^. For each value 
of ^, we compute E^g and E^g — E^^c for different values of the width parameter 
w. In Fig. 13-11 we plot E^g and E^q — i^vac as functions of C,, choosing the width at 
every point to minimize the energy (we do not allow w to be less than 1 so that we 
remain relatively insensitive to the Landau pole). For all points on the plot there is 
no spectral flow and so E^,^ = Ed + eiowcst + -f'vac in accordance with eq. ()3.23j) , where 
^lowest is the smallest positive bound-state energy in the Dirac spectrum. If i?vac is 
omitted, for < j9 < 0.6 we have configurations that have fermionic energies lower 
than the mass of the perturbative fermion (m^ = 1 in our units). These configurations 
indicate the existence of a local minimum on the E^,^ — E^^c surface which would be 
a soliton. The -Evac contribution, however, raises the energies of the configurations to 
above as shown in the figure, and the would-be solitons are destabilized. 

3.4.4 Paths over the Sphaleron 

The gauge fields introduce another mechanism for strongly binding a fermion state: 
there is a zero mode in the background of the sphaleron P7ll38j . Along an interpola- 
tion of the background fields from a configuration in C'-"^ to a configuration in C^'^~^^\ 
a fermion level leaves the positive continuum, crosses zero from above and finally 
enters the negative continuum. The lowering of the occupation energy, E^H, as we 
approach zero from above is offset by the rising effective energy {Ees = Ed + -Evac), 
so we must investigate whether the former can dominate the latter. We also use 
such interpolations to study the effects of a large Yukawa coupling on the sphaleron. 




(3.36) 
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Figure 3-1: Minimum fermionic effective energies (in units of mj), with as well as without 
E'vac contributions, along the interpolation in eq. p.36|) . 

We approximate the quantum-corrected sphaleron by minimizing the effective energy 
barrier between topologically inequivalent vacuum configurations, with respect to the 
variational parameters of our interpolations. We also investigate the possible emer- 
gence of new barriers in the one-fermion energy surface when the perturbative fermion 
becomes heavier than the quantum-corrected sphaleron. These last two phenomena 
affect the stability of the heavy fermion and in some models may be significant for 
baryogenesis. 

We make the following choices for the theory parameters: we fix the Yukawa 
coupling at / = 10, which is large enough that fermion effects are significant, but 
small enough to prevent our configurations from being affected by the Landau pole. 
Indeed we encounter no negative energy instabilities in our computations of E^^c 
for this coupling. We keep the Higgs mass fixed at v/\/2, which corresponds to 
m/j = ^7 ^ 0.07m f. We choose g to keep the quantum-corrected sphaleron energy 
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comparable to mj, since it is plausible that the mass of the sphaleron sets the scale 
for decoupling. If the sphaleron is too heavy compared to the fermion, the binding 
energy gained by the level crossing would be washed out by the effective energy 
of the sphaleron. For g = 6.5 our best approximation to the quantum-corrected 
sphaleron is approximately degenerate with the fermion. When the gauge coupling 
is given, the mass rriy^ of the gauge boson is determined from the renormalization 
constraint eq. ()3.17p : ^ 0.63m f for g = 6.5 and ^ 0.98m f for g = 10. These 
theory parameters are of course large deviations from the Standard Model parameters. 
We exaggerate them to see the effects of the heavy perturbative fermion. Another 
concern is that for large g, we should consider quantum fluctuations of the Higgs- 
gauge fields. However, we believe that anomaly cancellation drives the creation of a 
soliton, which would suggest that the fermion vacuum contains the essential physics, 
and our methods allow us to exactly compute this contribution to the energy for any 
Yukawa coupling. 

First we consider a linear interpolation between a winding-0 and a winding-1 
vacuum configuration. The interpolation parameter ^ goes from to 1: 

$ = ^,(1 _^)i + ^t,f/(i) ^ 
Wj = e^f/«9,[/«^ (3.37) 

In the spherical ansatz, U^^^ (x) = g{x) is specified by a single function, as in eq. ()3.27|1 . 
that we choose to be 

/(i)(r) = -27re~^/"', (3.38) 

where w characterizes the width of the configuration. We interpolate from ^ = (the 
trivial configuration with A^^cs = 0) to ,^ = 1/2 (a configuration with A^^cs = 1/2 in the 
presence of which the fermion has a zero mode) and vary w along the interpolation 
to minimize E^^^ . This gives an upper bound on the minimum E^^ as a function of 
A'cs- For A'cs = 1/2, this is an upper bound on the quantum-corrected sphaleron 
energy as well, because E^^^ = E^s in the presence of a fermion zero mode (the 
occupation energy is then 0). (Our numerical methods do not allow us to consider 
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Nqs exactly equal to 1/2, because the Higgs magnitude vanishes at r = and the 
second order Dirac equations develop a singularity as discussed in Appendix |X| But 
we can compute E'gg very close to Nqs = An exploration of Nqs between and 

1/2 is sufficient to map to all values of Nqs, since configurations with A^cs between 
1/2 and 1 are obtained by charge conjugation, and configurations with A^cs < or 
A'cs > 1 are large-gauge-equivalent to configurations with A'cs between and 1. 

We also consider an instanton-like configuration where the Euclidean time ^ = 
is the interpolation parameter (which varies from — oo to oo) between two topologi- 
cally inequivalent vacuum configurations: 



$ = y^Mr,Of/inst(x,0, (3.39) 



where 



is the canonical winding-1 map from (space-time infinity) to SU{2). Furthermore 
h is a function of the Euclidean space-time radius {y/r"^ + ^^) and goes from to 1 as 
this radius goes from to oo. 't Hooft's electroweak instanton j2] is constructed as 
a self-dual gauge field configuration in the topological charge one sector, and a Higgs 
field configuration that minimize the covariant kinetic term in the Lagrangian density. 
(In Sec. 12.2.41 this was constructed for the SU(2) theory with no Higgs fields.) This 
gives 



Hr,0= , (3.41) 

_|_ _|_ yjZ 



for any width w (the classical theory with no Higgs field is scale-invariant). We modify 
this radial function to exponentially approach its asymptotic value of 1, so that the 
potential in our Dirac equation falls off fast enough to have a well-defined scattering 
problem (as described in Appendix [X]). We choose 

/i(r,e) = 1-e- (3.42) 
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This choice does not minimize any part of the classical Euclidean action (in the 
topological charge one sector). Since we are interested in minimizing E^^ , which 
has fermion vacuum energy and occupation energy contributions in addition to the 
Higgs-gauge sector classical energy, we do not need our configurations to minimize 
the classical energy. In fact, as we describe later, the configurations that minimize 
are rather different from those that minimize Ed- 
In order to compute the Dirac spectrum for this background using the methods 
described in Appendix El we gauge transform to Wq = and limr^oo V = linir^oo = 
using the transformation function 

fir, = y_ de^:^^^,h{r^ e) - 27r (3.43) 

in eq. (|3.27|) . Finally, we go from ^ = — oo (the trivial configuration with Nq^ = 0) 
to = (a configuration with A^cs = 1/2 in the presence of which the fermion has a 
zero mode) and for each C, we choose the w that minimizes E^^g . 

In Fig. 13-21 we plot the minimum effective energies in both the zero-fermion sector 
(E^g) and the one- fermion sector (-Ees minimized within our variational ansatz 
for the linear interpolation, as functions of the Chern-Simons number, Nqs (see 
eqs. ()3.37I3.38|) ). As mentioned before, we fix the theory parameters at a Yukawa 
coupling of / = 10, a gauge coupling of g = 6.5 and a Higgs mass of m/j ^ 0.07m f. 
These parameters determine the mass of the gauge bosons to be ^ 0.63m/ from 
the renormalization constraint, eq. ()3.17|) . To isolate and highlight the contribution 
of the fermion vacuum energy, -Evac? we also plot the energies minimized with i^^ac 
subtracted. 

First consider the zero-fermion sector with the two curves E^^ and E'gg — E'vac- For 
all points on the plot there is no spectral fiow and so E^^ — E^^^c = Ec\. At A'cs = 0, 
both E^Q and E^i are minimized at the trivial vacuum configuration, eq. (j3.31|) with 
f(^)(^r) = 0. At A'cs = 1/2, E^^ is minimized at the quantum-corrected sphaleron 
while Ed is minimized at the classical sphaleron. Within our variational ansatz, 
we find the parameters that minimize E^,^ at Nqs = 1/2 are different from those 
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1.1 0.1 0.2 0.3 0.4 0.5 



^CS 



Figure 3-2: Minimum effective energies (in units of mj) along the linear path in eqs. (|.S.37[ 
I3.38|) . in both the zero-fermion and one-fermion sectors (with as well as without the -Evac 
contributions) . 

that minimize Ed- So our approximation to the quantum-corrected sphaleron is 
distinct from our approximation to the classical sphaleron. Moreover, the fermion 
vacuum energy correction to the sphaleron turns out to be rather large. Our classical 
sphaleron has an energy of 0.45mj (which agrees well with the numerical estimate 
of E = 1.52^^ in [IE]), while our quantum-corrected sphaleron has an energy of 
1.02m/. 

Next consider the E^^ and — -Evac plots in the one-fermion sector in Fig. 13-21 
Again, for all points on the plot there is no spectral flow and so E^^g = -E'ci+eiowcst+-E'vac 
in accordance with eq. ()3.23p . where eiowest is the smallest positive bound-state energy 
in the Dirac spectrum. Since the classical sphaleron has an energy much smaller than 
the perturbative fermion mass, one would expect that the perturbative fermion would 
have an unsuppressed decay mode over the sphaleron, as first pointed out by Rubakov 
in . The E^g — E^^^ curve indeed displays this decay path. The fermion vacuum 
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energy modifies things in two crucial ways. First, the fermion quantum corrections 
to the sphaleron raise its energy. (The theory parameters have been chosen so that 
the quantum-corrected sphaleron is approximately degenerate with the fermion, as 
mentioned before.) So the threshold mass is significantly increased. Second, in the 
plot of -Egg we observe that there is an energy barrier between the fundamental 
fermion and the quantum-corrected sphaleron. This indicates that even when the 
fermion becomes heavier than the sphaleron, there might exist a range of masses for 
which the decay continues to be exponentially suppressed (since it can proceed only 
via tunneling). 
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Figure 3-3: Minimum (in units of mj) for paths from a vacuum configuration to the 
sphaleron. Curves denoted 'Lin.' refer to the linear path in eqs. (|3.37l I3.38|) while those 
labeled 'Inst.' are associated with the instanton path, eqs. ()3. 3913.4013.42(1 . 

In Fig. 13-31 we restrict our attention to the one-fermion sector and consider the 
linear interpolation and the instanton-like interpolation. In addition \,o g = 6.5 we 
consider (7 = 10, which corresponds to 0.98m/. For each of these two gauge 

couplings, we plot the effective energy minimized within our ansatze as a function of 
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Ncs for the two interpolations. 

The two seemingly different interpolating configurations produce very similar min- 
imum E^g curves. Furthermore, when we enlarge the variational ansatze in our in- 
terpolations, we are unable to reduce the energies by any significant amount. Thus 
we speculate that the plotted curves may be close to tight upper bounds on the true 
minimum E^^ curve. This is the justification for considering only the linear interpo- 
lation in Fig. 13-21 and taking the evidence for the emergence of a new barrier and the 
significant energy change of the sphaleron seriously. 

Note that as the gauge coupling increases, lowering the energy of the quantum- 
corrected sphaleron, the barrier between the fundamental fermion and the sphaleron 
does not persist indefinitely. We observe that for g = 10, when is approximately 
1.3 times our quantum-corrected sphaleron, there is no barrier and the decay mode 
is finally unsuppressed. 

Just as in the case of the twisted-Higgs variational ansatz, in our ansatze of paths 
from a vacuum configuration to the sphaleron, we have not found a configuration 
with the associated fermion energy lower than both the perturbative fermion and the 
quantum-corrected sphaleron. Thus, we find no evidence for the existence of fermionic 
solitons in the low energy spectrum of the Standard Model within our ansatze. 

3.5 Conclusions and Discussion 

We have explored quantum effects of a (hypothetical) heavy fermion on a chiral gauge 
theory (the electroweak theory, ignoring the hypercharge gauge fields) for certain 
parameters. The quantum-corrected sphaleron is heavier than the classical sphaleron 
by an energy of the order of the perturbative fermion mass. This higher barrier 
suppresses fermion number violating processes. We also observe the creation of an 
energy barrier between the perturbative fermion and the sphaleron, so that a fermion 
with energy slightly above the sphaleron can still only decay through tunneling. This 
new barrier exists only for an intermediate range of perturbative fermion masses, and 
a heavy enough perturbative fermion is indeed unstable. We do not, however, find 
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evidence for the existence of a soliton in the spectrum of the theory. The fermion 
vacuum contribution seems to destabihze any would-be sohtons. 

In hindsight it is not too surprising that we do not find a quantum sohton within 
the spherical ansatz. As mentioned in Sec. 13.11 the existence of such solitons would 
maintain anomaly cancellation when fermions are decoupled from the electroweak 
theory. Without the hypercharge gauge field, the only anomaly is Witten's global 
anomaly j21j due to topologically non-trivial maps from 5*4 to SU{2). However, in 
the spherical ansatz, the theory reduces to a U{1) theory in which Il3{SU{2)) persists 
as ni(f/(l)) (and so we have topologically inequivalent vacua and the weak sphaleron) 
but there is no remnant of Il4{{SU{2)). So, the quantum soliton that could resolve 
the decoupling puzzle probably lies outside the ansatz. 

One route beyond the spherical ansatz is the use of non-trivial Il4{SU{2)) to 
construct non-contractible loops and corresponding novel classical sphalerons (see 
Sec. I2.2.6p . We expect fermion zero- modes in these sphaleron backgrounds and thus 
the required tight-binding-mechanism exists. Moreover, since the non-trivial topology 
that gives rise to Witten's anomaly is built into the construction of such configura- 
tions, they are promising candidates for objects that maintain anomaly cancellation 
in the decoupled theory. However, an evaluation of the stability of such a state would 
require a calculation of the one-loop effective energy, which is more difficult in the 
absence of spherical symmetry. We have not pursued this direction any further. 

Another interesting non-spherical configuration is the W-string solution discussed 
in Sec. 12.2.41 Our motivation goes beyond fermion decoupling. Stable electroweak 
strings are a crucial ingredient in a scenario for electroweak baryogenesis even without 
a first order electroweak phase transitional]. It is likely that the classically unstable 
W-string could be quantum-stabilized by having it carry fermion number. The string 
configuration gives rise to fermion zero modes and it is possible that that several 
massless quarks and leptons are trapped along the string in a stable manner. We 
discuss this in greater detail in the next chapter. 
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Chapter 4 



Electroweak Strings 



In this chapter we begin by reviewing the Higgs-gauge sector of the electroweak the- 
ory, after reintroducing the hypercharge gauge fields that we have ignored so far. We 
then describe the electroweak string solutions in the theory by constructing them as 
embedded Nielsen- Olesen vortices. We explain the significance of these objects and 
how they may be stabilized by trapping heavy quarks. An investigation of this sta- 
bilizing effect of quarks requires us to compute fermion vacuum energies in the back- 
grounds of strings. Our computational method involves scattering data of fermions. 
However, the string backgrounds generate long-range potentials leading to difficulties 
with scattering theory. In order to isolate and examine these issues, we study the sim- 
pler problem of magnetic flux tubes in QED, as a stepping stone to the electroweak 
strings analysis. 



4.1 The Higgs-gauge Sector 

In Sec. 12.11 we described a simplified version of the bosonic sector of the electroweak 
theory in which we ignored the hypercharge gauge fields, B^. Now we re-introduce 
them and consider the full SU{2) x U{1) gauged Higgs theory. It is defined by the 
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action 



d X 



(4.1) 



The non-Abelian field strength tensor, W^^, is defined in Sec. 12.11 The Abehan field 
strength tensor is 

(4.2) 



B^u = d^By — dyBi^ . 



The Higgs doublet is 



(4.3) 



and the covariant derivative is 



D.ct) = {d^ - igW,^ - ig%B,) 



(4.4) 



where the Higgs hypercharge is = |. Under a gauge transformation V = Vi x V2, 
where Vi = exp{ip{x)Y^) E U{1) and V2 = exp(m''(x)rV2) G SU{2), 



W„ 



^1^20, 

B^ + ^d^P, 
9 

V2 (w^ + -d,] 



(4.5) 



After spontaneous symmetry breaking, when the Higgs acquires a vev (($) = 
vl), the gauge fields eat the three Goldstone bosons and acquire masses. The mass 
eigenstates of the gauge fields (in unitary gauge) are 



A, 



smO^Wj^ + cose^Bf,, 



cos6'^lV^ - sin 9i,B^ , 



W, 



± 



72 ^ ^ ^ ' 



(4.6) 



58 



where the mixing angle is defined as 



tan d,, 



9 



(4.7) 



The photon is massless (due to the unbroken f/(l) subgroup of electromagnetism), 



the neutral boson has mass mz = '^yi.O'Y + 9^'^i the two charged W bosons 
are degenerate with mass mw — '^9'^- The neutral Higgs field 0o has mass 
niH — 2v-\/A. In terms of the mass eigenstates, the action becomes 



Sh[^.Z,,W^,A,] = jd'x 



+ (l>'^0)^l>^0-a(</)V-^')' 



(4.8) 



The field strength tensors are Abehan, the covariant derivative is 



(d, - t^Jg^ + g'^QzZ^ - icQaA, - i-^ (w^T+ + W^T-)^ <t> , (4.9) 



where 



e 

Qa 
Qz 



99 



T^ + Y, 



cos\e^)T^ - sm\e^,)Y = - sin2(^^)g^ , 



(4.10) 



We use the experimental values for the mixing angle and the fine structure constant 
{sin^ 9w — 0.23, Q! — 1/137) to determine 



9 = 



9 



sin 9,, 



cos 6„ 



0.63, 



0.31 



(4.11) 
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The observed Z mass of 90 GeV fixes the vev at f ~ 177 GeV. Finally, choosing the 
Higgs mass to be 115 GeV we get the self-coupling A ~ 0.11. 

4.2 The String Solutions 

Now we consider static configurations that are trivial in one direction, say z, and have 
finite energy per unit length in a localized region of the x — y plane. Since the region in 
space in which the energy is localized extends in the z direction, such configurations 
are called strings. They are also called cosmic strings to distinguish them from 
fundamental strings in string theory. Finite energy requires the configuration to 
be pure gauge at planar infinity. Let (p, ({>) denote the radial and angular coordinates 
on the plane. If Vi, V2 are elements of t/(l), SU{2) respectively, then the asymptotic 
configuration at p ^ 00 is defined by a map Vi{ip)V2{ip) from 5*1 (the planar boundary 
spanned by Lp) to the gauge group: 

Wlr^ = '-V2d,V^. (4.12) 

If the gauge group had been U{1), then these maps would be characterized by wind- 
ing numbers {Ili{U{l)) = Z), and we would obtain topologically stable Nielsen- 
Olesen vortex solutions [Hj. (They are called vortices because the gauge fields are 
azimuthally directed, as we shall see later.) However, the gauge group is SU (2) x f/(l). 
But we can still embed Nielsen-Olesen solutions in this theory by considering maps 
from the planar boundary to various U{1) subgroups of the gauge group that break 
completely when the electroweak symmetry breaks. These solutions are unstable. 
This is readily seen if the string solutions are constructed using non-contractible 
loops of configurations, as was done in the SU{2) theory in Sec. 12.2.41 In that case 
it was shown that the strings are sphalerons with two directions of instability. For a 
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comprehensive review of electroweak strings, see |lUj . 



4.2.1 The Z-string 

First consider the U{1) subgroup generated by the charge Qz associated with the 
gauge bosons, as defined in eq. I4.1U1 The maps that specify the asymptotic pure-gauge 
configurations are 

V{^) = e-2™^Q^ , (4.13) 

where n is an integer (otherwise the Higgs fields are not single-valued under (f — 
(f + 2tt). Then, the winding n Z-string solution is 

00 = /^^(p)^e^"^ 
2n 

Vg + 9 p 

with all the other fields vanishing. The unit vector in the azimuthal direction is 

(4.15) 




The radial functions fn-, fc vanish at the origin (for smooth fields) and approach 1 at 
infinity (for finite energy per unit length). They satisfy the Nielsen- Olesen differential 
equations 



1 r?^ 
P P 

= fi--f'G + W + 9>'fHa-fG), (4.16) 
p ^ 

and can be determined numerically. 

The Z-string is a vortex configuration carrying Z magnetic fiux of 4?™/-*/ + g''^. 
The magnetic pressure tends to spread out the fiux. However, the Higgs condensate 
is suppressed in the region of the Z magnetic field, and this tends to compress the 
fiux. The two competing effects result in an equilibrium size for the Z-string. How- 
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ever, had we considered the unbroken electromagnetic U{1) subgroup generated by 
Qa: the Higgs condensate would be unaffected in regions of magnetic field, and the 
configuration would be driven to a singularity. (The energy could be made closer to 
zero by increasing the size of the area of the fiux and independently decreasing the 
area in which the Higgs condensate is suppressed.) 

4.2.2 The W-st rings 

Consider the 1-parameter family of U{1) subgroups generated by 

r^ = cos(OTi + sin(Or', (4.17) 

where ^ is between and tt. The maps that specify the asymptotic pure-gauge 
configurations are 

V^{ip) = e^"'^^« (4.18) 
where n is an integer. This gives a 1-parameter family of winding n W-string solutions: 

/ ie-'^sm{nip) \ 

(P = Jh[p)v 

\ cos{n(p) I 

= _cos(0/g(p)— 
9 

2n 

W = -sin(0/G(p)y<^. (4.19) 

The radial functions satisfy the Nielsen- Olesen differential equations. 

In analogy with the Z-string, W-strings are vortex configurations carrying W^"^ 
magnetic flux of ^-nnjg. Their size is stabilized by the competition between the 
magnetic pressure and the potential energy. 
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4.3 Fermions on Strings 

Although the electroweak strings are classically unstable, they may be stabilized by 
their interactions with quarks and leptons. Since the Higgs condensate is suppressed 
in the core of the string, the fermions become massless on the string, and could 
resist its decay. Of course, as discussed in Chap. El to be consistent to order Ti, the 
fermion vacuum energy contribution must also be included in the stability analysis. 
The minimum total energy associated with Nf fermions (with perturbative mass nif) 
trapped on the string is 



with the different energies above defined in Table IXTl It is conceivable that for some 
value of A'^^, 



leading to stable multi-quark objects with possibly rich phenomenology. 

Electroweak strings, if stabilized by the fermion binding mechanism, could have 
profound cosmological consequences. A gas of strings would have negative pressure 
associated with it (the strings tend to contract in length) and could contribute to 
the dark energy that is required to explain the recently observed cosmic acceleration 
[121 US] • Also, as pointed out by Nambu jlOl , Z-strings are expected to terminate in 
monopole-antimonopole pairs, which could give rise to a primordial magnetic field. 
Furthermore, a network of stable strings at the electroweak phase transition would 
provide a scenario for electroweak baryogenesis without requiring a first-order phase 
transition ^T]. As the strings contract, they provide out-of-equilibrium regions, and 
the core of the string has copious baryon number violation due to the suppressed Higgs 
condensate. This is an alternative to the usual idea of bubble-nucleation baryogenesis 
which requires a first order phase transition to go out of thermal equilibrium. 

Our calculational method for the fermion vacuum energy uses scattering data of 
fermions. This raises some technical issues because the string configurations give long- 



(4.20) 
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range potentials with subtleties associated with scattering theory. So, as a stepping 
stone to coupling fermions to electroweak strings, we have considered the simplest 
example of a vortex configurations: magnetic fiux tubes in quantum electrodynamics 
(QED). This allows us to isolate and examine the problem of scattering fermions off 
vortices. We discuss this in detail next. 

4.4 Flux Tubes in QED 

We investigate quantum corrections to the energy of a magnetic fiux tube in QED, 
up to one-loop order. We consider space-time dimensions of D=3 and D=4. In D=3 
the configuration is a vortex (rather than a tube). In D=4 we consider the energy per 
unit length of the fiux tube. Our primary motivation is that this analysis will shed 
light on vortices (strings) in more complicated field theories, especially the Z-string 
in the standard electroweak theory. 

We compare the energies in three and four space-time dimensions to see if the sim- 
pler case of D=3 (with no divergences) bears any resemblance to the D=4 calculation. 
The classical energy is of course the same in the two cases. The quantum correction 
to the energy could possibly be very different, especially because of the different di- 
vergence structure in the two cases. In D=4, the bare one-loop energy is divergent 
and only after renormalization conditions are imposed on the calculation do we get 
a finite result. In D=3, in contrast, the bare energy is finite. But, the comparison 
between the two cases is sensible only when we use the same renormalization condi- 
tions. This requires a finite counterterm in the D=3 case to keep the renormalized 
photon field and electric charge fixed. Without this finite renormalization, the D=3 
and D=4 energies are qualitatively different. But, after proper renormalization, we 
find that the energies are closely related. This suggests that if quantum corrections 
stabilize the Z-string in two spatial dimensions, then that might be indicative of a 
stable Z-string in three spatial dimensions. 

We also use this problem to get a handle on several technical issues associated with 
the computation of the one-loop energy of a fiux tube. As described in Sec. 14.4.21 an 
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efficient way to compute the energy is to use scattering data of fermions in the back- 
ground of the flux tube. However, vortex configurations give long range potentials, 
which violate standard assumptions in scattering theory |41j. This leads to inconsis- 
tent and difficult to interpret results, some of which are discussed in Sec. 14.4.31 We 
show that these puzzles are due to the fact that an isolated flux tube is unphysical, 
and once a region of return flux is included to make the net flux 0, the scattering is 
well-defined and the puzzles disappear. In the limit of the return flux being infinitely 
spread out, the energy has a well defined limit, which is associated with the localized 
flux tube. 

4.4.1 The Theory 

We consider QED in space-time dimensions D=3 and D=4, with a four-component 
fermion field, ip, in both cases. The Lagrangian density is 

= -\f^uF^'' + e4- m)ij + , (4.22) 

where the indices run from to D-1 and is the counterterm Lagrangian in 

D space-time dimensions. In D=4, this describes electromagnet ism with a single 
fermion flavor of mass m and charge e. In D=3, we have parity-invariant electro- 
magnetism with two flavors of fermions (of equal mass m and equal charge e). The 
four-component fermion contains two two-component fermions. This can be thought 
of as a planar version of the D=4 theory in which all particles (including virtual ones) 
are constrained to lie on the x-y plane. 

We are interested in static magnetic flux tubes. These are localized, cylindrically 
symmetric magnetic fields (pointing in the z direction in D=4), with a net flux F 
through the x-y plane. For example, a Gaussian flux tube of width w is 

Bg{p) = Bo{Q)e-P''-' , (4.23) 
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where — -\-y^. Its flux is 



Fq = 'kw'^Bg{{)) . (4.24) 

We often find it convenient to specify the fiux in units of 27r/e and define 

F = JT— . (4.25) 
e 

In the radial gauge, vortex configurations of gauge fields give fiux tubes: 

Ao = 0, 

A = ^^np)0 (4.26) 

where f{p) goes from to 1, and 

(4.27) 

For small p, /(p) oc p^, so that S is not singular. 
4.4.2 The Energy 

We want to compute the energy, E^^\ (per unit z- length in D=4) up to one- loop order 
of magnetic fiux tubes. For static configurations, this is minus the one-loop effective 
action per unit time (per unit z- length in D=4). Since the theory is Abelian, the 
photons do not self-interact and the one-loop effective action is obtained by integrating 
out the fermion field. The photon fiuctuations start contributing only at two loops 
and higher, and we ignore their contributions. Then, 

= E,, + EiZ^ (4.28) 
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where the classical energy (in both D=3 and D=4) is 

Eci^lJ d^xB-" , (4.29) 

and the renormalized fermion vacuum energy (the renormalized one-loop energy) is 

EiZ^ = ^\im^ ^ [In Bet{tp+ eyl- m) - In Det(z^- m)] + . (4.30) 

The terms in {} above are present only in D—A where the energy per unit z-length is 
required. The counterterm energy, , and the functional determinant contribution 
to E^^J are discussed in the following subsubsections. Note that the energy is defined 
relative to the absence of electromagnetic fields. 

Renormalization 

The fermion fluctuations renormalize the bare photon fleld and fermion charge, 

e = (l + C(^))^/'eba.e, (4.31) 
giving the counterterm energy 

E'g^ ^C^"")]^! d^xB\ (4.32) 

In the absence of photon fluctuations, the bare fermion fleld and mass do not get 
renormalized. In D = 4 the counterterm is divergent and combines with the diver- 
gent functional determinant to give a flnite, renormalized one-loop energy. In D — 
however, the renormalization of the bare photon fleld and charge is flnite resulting in a 
flnite counterterm. Nevertheless, the counterterm must be taken into account, other- 
wise the photon fleld is not properly normalized and the charge is the (unobscrvable) 
bare charge. Renormalization is the process of flxing the unobscrvable bare parame- 
ters in the theory by imposing physical conditions. In D=4 this process also removes 
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the ultraviolet divergences (by allowing the bare parameters to be divergent), but the 
lack of divergences in D=3 does not mean that renormalization is not required (as is 
often stated in the literature). In fact, one of our aims is to compare the energies in 
D=3 and D=4, and this is sensible only when we impose the same renormalization 
conditions in the two cases. 

We choose the on-shell renormalization condition (residue of the = pole of 
the photon propagator is 1) to fix the counterterm coefficient: 

C(^) - r (2 - f ) 

^ " 3(47r)^/2 m^-D ' ^^'^^^ 
In D = 3, it is the finite quantity 

C^^) = -^, (4.34) 
and in D = 4 it is the divergent quantity 

Note that we have dimensionally regulated the divergent counterterm coefficient in 
D=4 by introducing the regulator e = A — D. 

Functional Determinant 

We use the phase shift approach, described in Sec. 13.2.21 to exactly compute the 
functional determinant that contributes to the one-loop energy as in eq. I4.3U1 Recall 
that the basic idea is that the fermion vacuum energy is given by the sum over the shift 
in the zero-point energies of the fermion modes due to the presence of the magnetic 
background. In the continuum this is a sum over bound state energies and an integral 
over the continuum energies weighted by the change in the density of states. Since 
the vortex configuration that characterizes the fiux tube is cylindrically symmetric, 
the scattering matrix of the fermion may be decomposed into partial waves labeled 
by the z-component of the total angular momentum, M. Then, in D=3, the change 
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in the density of states is given by the phase shifts of the fermion scattering wave 
functions in the flux tube background: 

where k is the magnitude of the momentum. The phase shifts above are defined as 
the sum over both signs of the energy. In D=4, there is an additional factor in the 
change in the density of states corresponding to the trivial z dimension, and we use 
the interface formalism [321 to straightforwardly account for this. We subtract the 
first two terms in the Born series expansion of the change in the density of states and 
add the corresponding energy back in the form of the two-point Feynman diagram 
(a fermion loop with two insertions of the background potential). In D=4, the Born 
subtraction renders the integral over the continuum energies convergent and isolates 
the ultraviolet divergences in the Feynman diagram, which when combined with the 
counterterm gives a properly renormalized finite result. In D=3 we did not have to 
Born subtract, but we do so anyway to maintain a close similarity with D=4. In 
Appendix El we provide the details on how we compute the phase shifts and the born 
series. 

Our final expression for the renormalized fermion vacuum energy is 

Ei^J = 4? + EiS + Ef . (4.37) 

The unrenormalized Feynman diagram energy (unrenormalized two-point energy) is 
87rJ^2 .oo f r , df{p) ^ . y f\ x{l-x)pT{2- D/2) 



'''^ = (4^ Jo ^nio '^P^-^'^^'Pn Jo 



(47r)^/2 Jo \Jo dp ''^"^'J Jo (m2 + p^il - x)^-^/^ ' 

(4.38) 

which combines with the counterterm energy to give the renormalized Feynman di- 
agram energy (in D=4 this cancels the divergence). The phase shifts energy, which 
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corresponds to the contributions from the three-point and higher functions, is 



m 



M 



^1'^ = I dkkXn ^'' ' Y.'^M{k), (4.40) 

27r^ Jo m ^ 



2 



where 



5M{k) = 6M{k) - 4/ (^) - ^M{k) , (4.41) 



and denotes the i^^ term in the Born series expansion of the phase shift. 



4.4.3 A Few Puzzles 

In the previous subsections we have described how to compute properly renormahzed 
energies of magnetic flux tubes, to one-loop. However, before proceeding with the 
calculation, there are several subtleties that need to be addressed. The main issue 
is that even for exponentially localized magnetic fields, the photon field falls off too 
slowly (like 1/p in the vortex configuration) for the standard theorems of scattering 
theory to hold. Since our computational method relies on scattering phase shifts, 
we need to resolve this. In addition to scattering theory problems, there seem to be 
rather unexpected and physically senseless results associated with flux tubes. In this 
subsection we enumerate some of these. 

As described in Appendix El we use the second order Dirac equations in two radial 
functions, fi'i(p), fi'2(p), to calculate phase shifts in the background of a flux tube. The 
asymptotic form of these equations is 

^ ,'iJj^(t'-^Ji±i±Il),.. (4.42, 

So, instead of the free equations, we have the equations in the presence of an ideal 
flux tube (with profile function /(p) = 1 and singularities at p = 0). The centrifugal 
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barrier is shifted by the flux and the regular solutions are Bessel functions: 



9l{p) = J\M^\+T\{kp),g2{p) = J|A/+i+^|(M • (4-43) 

From the asymptotic forms of these [kp » 1), we find a phase shift relative to the 
trivial configuration (A^ = 0) by the amount 

SM,±,tee{k) = ^ (|M - i| + |M + i| - |M - i + ^1 - |M + i + ^1) , (4.44) 

where the second index labels the sign of the energy oo ^Sj- The second order equations 
allow us to compute phase shifts relative to the ideal flux tube, 5m,±, ideal, to which 
'^M,±,free must bc added to obtain the phase shift relative to the trivial configuration: 

Sm,± = SM,±,idea.\ + ^M.i.free • (4.45) 

However, for k = 0, there is no p large enough where the asymptotic form of the 
Bessel functions holds, and it is unclear what 5M,±,free(0) (and hence 5a/,±(0)) should 
be. 

Now let's restrict our attention to momenta strictly greater than 0. In a delta- 
shell background of magnetic flux at p = po, the phase shifts can be calculated 
analytically We consider more general localized flux tubes and compute phase 
shifts numerically. We agree with the analytical calculation on the expressions for 
the phase shifts at large and small momenta. These expressions depend only on the 
net flux through the plane and not on the particular profile of the magnetic field. In 
the Jz = M channel, at large momentum, 

lim 5M,±(fc) = 0. (4.46) 

For small momentum, 

f 7r(|M + i| - |M + i + ^|) : M> i 
i 7r(|M-i|-|M-i+^|) : M < -i 
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Armed with the hmiting values of the phase shifts at large and small momenta, and 
using the fact that the momentum derivative of the phase shifts gives the density of 
states, we can ask how many states leave the continuum as the flux tube is turned on. 
This should correspond to the number of bound states. For concreteness, consider 
T = 1.2. In this background there are two bound states: one with = 1/2 and 
u = —m, and another with = —1/2 and uj = m. However, the number of states 
that have left the continuum in channel = M is 



So, in no channel is there an agreement between the number of states that leave the 
continuum and the number of bound states and we have a puzzling discrepancy in 
the number of states. 

Apart from the formal scattering theory problems discussed above, there seem 
to be physical puzzles involving a flux tube. These show up in D=3 with a single 
two-component fermion. In this case, the net fermion charge in the background of 
the flux tube is as is well-known (see jlHI for example). This implies that as 

the flux tube is turned on, charge conservation is violated. Even if we ignore this 
problem of how to generate the flux tube without violating charge conservation and 
assume that the flux tube is given as an initial condition, we run into a problem. As 
the flux tube is spread out so that the magnetic field approaches everywhere (but 
the flux is held fixed), the total energy of the flux tube approaches (as shown in 
Sec. I4.4.5|) . This means that we have charged, massless objects in the theory and the 
fundamental fermion cannot be stable. 

4.4.4 Embedding 

All the subtleties and puzzles discussed in Sec. 14.4.31 could perhaps have resolutions. 
However, we now argue that it is impossible to create a configuration carrying non-zero 



1.6 



M 



1 

2 




M<-1 



(4.48) 
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net flux starting from nothing. We demonstrate how the flux tube can be embedded 
in a physical no-net-flux configuration in which the localized flux tube has a spread 
out region of return flux. Now the scattering potential falls into the class of potentials 
well understood in scattering theory and all the associated puzzles disappear. As the 
return flux is spread out, we show that the energy of the embedded configuration 
approaches a well-defined limit which we interpret to be the energy of the flux tube. 
This is exactly the energy obtained by ignoring all the problems with the phase shifts 
in the background of a flux tube and using them to compute the energy. 
In QED3_|_i we have the Bianchi identity 

e^^'^'^a^jF^.^O. (4.49) 

This is a mathematical identity and not an equation of motion from the action prin- 
ciple. In terms of electric and magnetic fields it gives us the two Maxwell equations: 

V-B = 0, 
— = -Vxi. (4.50) 

The fact that the magnetic field is divergence-less requires all magnetic flux lines to 
be closed. So, there can be no net flux through the x — y plane. In QED2+1 , only 
the second of the two equations remains, in the form 

— ^-d,Ey + dyE,. (4.51) 

Integrating both sides over the area of the x-y plane, we get: 

F(t) = F(i = 0) - / E- dt, (4.52) 

J c 

where C denotes the circle at p ^ 00. So, if we start with net flux zero, we can create 
a non-zero net flux only if we have an electric field at spatial infinity (with non-zero 
curl), which is physically impossible. Thus, both in D=3 and in D—4, every flux tube 
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that is set up locally must have a compensating return flux to make the net flux zero. 
This return flux is not a crutch to enable the use of scattering theory to compute 
energies, but rather a physical requirement. 

The presence of the return flux not only removes the subtleties associated with 
scattering theory, but also eliminates the physical puzzles discussed in the single 
fermion D=3 theory at the end of Sec. 14.4.31 If we can't create a net flux, we can't 
create a net charge either. We can only polarize the Dirac sea and create local regions 
of charge with net charge zero. Furthermore, by requiring the net flux to be zero, we 
have eliminated the presence of massless charged objects in the theory. 

Now we demonstrate that as the return flux is spread out so that the associated 
magnetic fleld, B^, locally approaches 0, the energy approaches a well-deflned limit. 
This can be understood analytically by expressing the total energy as a series in 
inverse powers of the width of the magnetic fleld. As shown in Sec. 14.4.51 for a flxed 
flux, as the magnetic fleld is spread out, the classical energy goes to like l/w"^ and 
the renormalized one-loop energy goes like l/w"^ (in both D=3 and D=4). Thus, in 
the limit of spread-out return flux, the energy should limit to that associated with the 
local flux tube alone. We show that this value corresponds to computing the energy 
using the phase shifts in the unphysical problem with no return flux. 

Consider the Gaussian flux tube, Bq{p), deflned in eq. 14.231 with flux F. Choose 
the return flux to be centered at a ring of radius pn (with a width proportional to 
Pr), e.g. 

Br{p) = 7 4.53 

Trp% (l -I- 256 (pVpI - 1)^ j (7r/2 + arctan 16) 

The local flux tube is superimposed with the return flux to get the no-net-flux em- 
bedding 

Boip) = Bg{p) + Br{p) . (4.54) 

The limit of spread out return flux corresponds to p/j ^ oo. The total energy consists 
of three parts: the classical energy, the renormalized Feynman diagram energy and 
the phase shifts energy. We consider these different energy contributions for Bq as a 
function of p/j. We scale all dimensionful quantities in units of the fermion mass, m 
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Figure 4-1: Renormalized Feynman Diagram energies (in D=3 and D=4) as a func- 
tion of the return flux radius. The sohd hues correspond to the energies without the 
return flux. 

and choose JF = 4.8 and w = 1. Using the expression for the classical energy, it is 
trivial to see that 

lim E,i[Bo] = E,i[Bg] . (4.55) 

In Fig. 14-11 we show the approach of the renormalized Feynman diagram energy, 
+ Eqt, in the no-net-flux configuration to that in the flux tube configuration, 
as the return flux is spread out. In both D=3 and D=4 there is good convergence 
for pn > 10. In Fig. 14-21 we show the integrand of the phase shift part of the energy, 
in both the embedded configuration (for two fixed values of pn) as well as in the flux 
tube configuration, for D=3. Note that the integrands disagree only for small values 
of the momentum (for which the scattering fermion is sensitive to the presence of 
the spread out return flux). The integrand in the embedded configuration oscillates 
around the integrand in the flux tube configuration, with an amplitude that decreases 
as k increases. The region of disagreement gets pushed to smaller values of k as pr 
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Figure 4-2: Integrand of the phase shift contribution to the energy in D=3 for pn = Q 
(left panel) and pr = 26 (right panel). The solid lines correspond to the isolated flux 
tube problem and the dashed lines to the embedded problem. 

increases. In the spread-out limit, we find that the phase shift contribution to the 
energy is identical in the embedded problem and the flux tube problem. We have 
verified that the same results hold in D=4 and do not plot the corresponding results. 

Thus, when the flux tube is embedded in a physical no-net-fiux configuration and 
the return flux is spread-out, we find that the energy in the embedded problem ap- 
proaches the energy in the isolated flux tube problem. This allows us to use the phase 
shifts in the background of the isolated flux tube to compute the energy, with the 
implicit understanding that there is a spread-out return flux present, which eliminates 
all the puzzles and does not contribute to the energy. 

4.4.5 Results 

Derivative Expansion 

For a fixed magnetic field at the origin, we vary the width of the flux tube (thereby 
varying the flux as well). As the width increases, the configuration approaches a 
non- vanishing constant magnetic field, and the first few terms in the expansion of the 
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energy in derivatives of the magnetic field should be a good approximation. In this 
subsection we investigate the convergence of the derivative expansion for the Gaussian 
flux tube of eq. KT^ 

In D=3, the derivative expansion of the unrenormalized one-loop energy is given 
by jiniiTl 

-^vac -^CT ~ -^vac.O ^ -^vac,2 5 l^^.dUJ 

to lowest non-trivial order in the derivative, where 

E'l, = ; f d'x\t{eB f\i^eBr'' f fe-"'/l'^l.-'/^ '''<7f ^' . (4.57) 

Note that the above expressions are twice those for a two-component fermion. Since 
our Gaussian flux tube is a function of p/w, we can rescale the x-y coordinates in 
units of w and obtain -Eyac.o ^ ""^^ -E'vac,2 w^. On dimensional grounds, it is 
clear that the counterterm contribution to the energy is proportional to w"^ and for 
the Gaussian flux tube. 

We add the counterterm to get the renormalized derivative expansion 

EfX = Egl + 42,0 + Ei% , (4.59) 

where the first two terms are proportional to w'^, the last term is independent of w, 
and all omitted terms go to for large w. 

In D=4, the lowest non-trivial order derivative expansion of the renormalized 
one-loop energy is [IE] 

Eit^c = Eito + 4£,2 , (4.60) 

where 

\eB\^ /-oo __2/,„o, _o / , , . 1 s 



Eito = Jd\^-^fdse—'/\^^\s~'[coth{s) 



s 3J ' 
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Figure 4-3: Renormalized one- loop energies in D=3 (left panel) and D=4 (right 
panel), for fixed values of the magnetic field at the origin, as a function of the width 
of the Gaussian flux tube. The lines correspond to the derivative expansion to lowest 
non-trivial order in the derivative. From top to bottom, eBc^O) = 1.1, 2, 2.5, 3, 4, 5. 



(4) 
vac, 2 



J d^x\V{eB)\'^\327r^eB\-^ dse-'""'/^^^^ (l - 4 coth^ s + 3 cothS+ 
;i -coth^s) ) . (4.61) 



3 coth s 



Again we rescale the coordinates in units of w and obtain -Ejac,o ^'^d -Evac,2 

In Fig. I4-3| we compare the exact -E^^^ with the renormalized derivative expansion 
approximation, for different fixed values of 6-80(0) and find excellent agreement for 
widths larger than 1. There appears to be no qualitative difference between the 
energies in D=3 and D=4. This is in contrast to the claim in [IHl that in D=4, 
renormalization effects cause the convergence to the derivative expansion result as a 
function of the width to be slower. Had we not included the counterterm in D=3 
(thereby having a different renormalization condition from the one in D=4), the 
corresponding energies would have been positive. 
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Fixed Flux 

Now we consider flux tubes with fixed flux and varying widths. We show that in the 
limit of the width going to infinity, the energy of the flux tube goes to zero. 
First consider the classical energy: 



e^w"^ Jo X \ dx^ 



g2^2 



for a Gaussian flux tube , (4-62) 

(4.63) 



where x = p/w. For large widths, this goes to like In the large width limit, the 

magnetic field becomes weak and so the dominant contribution to the one-loop energy 
comes from the two-point function. In D=3, the unrenormalized two-point energy 

can be expressed as a series in The leading term proportional to l/w"^ turns 

out to be exactly equal to minus the counterterm energy and only the sub-leading 
term remains in the renormalized two-point energy. For a Gaussian flux tube. 



w 



(4.65) 



So, for large widths, renormalization not only changes the sign of the one-loop energy, 
but also the rate at which is approached. The renormalized results in D=3 are 
similar to the results in D—A: 

41^ + 4'^ oc (4.66) 

to leading order in l/w"^. For a Gaussian flux tube the proportionality constant is 
l/(307r). Note that in the w ^ oo limit, the total energy in both dimensionalities 
is 0. This is why in our embedding of the flux tube in a a no- net-flux configuration. 
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Figure 4-4: Renormalized fermion vacuum energy (normalized in units of JF^) as a 
function of the width, for various fixed values of the flux T in the Gaussian flux tube 
{T = 2.5,4.5,6.5,8.5 from top to bottom). The left panel is for D=3 and the right 
panel for D=4. 

the energy approached a well-defined limit as the return flux was spread out (see 
Sec. 11331) • 

In Fig. 14-41 we plot the exact one-loop energies for various values of the flux as 
a function of the width. We normalize the energies in units of JF^ so that all differ- 
ences between the different fluxes are due to three-point and higher contributions. 
We compare the energies with the leading order two-point energies. We find good 
agreement for large widths. As the flux increases, we need to go to larger widths to 
get a weak magnetic field everywhere, and so the convergence to the leading order 
two-point energy occurs at larger widths. 



4.5 Conclusions 

We have argued that electroweak strings could trap quarks and give rise to novel 
objects with rich phenomenology. These could be very significant for electroweak 
baryogenesis and could also be dark energy candidates. In order to gear up for the 
electroweak strings stability analysis, we studied the simpler problem of fermions 



80 



in the background of magnetic flux tubes in QED. We found that all puzzles and 
problems associated with the long-range nature of the potential disappear when we 
consider a region of return flux (which is the physical scenario). This suggests that the 
electroweak strings calculation could be more tractable if we embed the background 
into one with no net flux. Also, we find that when the energies are properly renormal- 
ized, there are no qualitative differences between two and three spatial dimensions. 
So, it may be sensible to consider electroweak vortices in two spatial dimensions, and 
look for a stable fermionic soliton (in analogy with the quantum solitons discussed in 
Chap. El). If we find such an object, then it may suggest that a stable object exists 
in the physical case of three spatial dimensions. We are currently investigating this. 



81 



Chapter 5 



Conclusions and Summary 



We have used the non-trivial topology of maps from Ss and 5*4 into the gauge group 
SU (2) of the weak theory to construct non-contractible loops of Higgs-gauge config- 
urations in Euclidean spacetime. A winding 1 map from 5*3 to SU{2) allows us to 
construct the well-known weak instanton, weak sphaleron and the family of W-strings. 
When the winding is chosen to be an integer n, the method suggests the existence 
of muUi-instantons (with topological charge n), multi-sphalerons (barriers between 
winding vacuum and winding n vacuum) and vorticity n W-strings. However, these 
do not produce qualitatively different physical effects from those of the n — 1 solu- 
tions. For example, the n-instanton would mediate the production of n fermions via 
quantum tunneling through the sphaleron barrier. But {n times) repeated instances 
of the 1-instanton already allows for this process, albeit with a possibly different 
probability amplitude. A topologically non-trivial map from 5*4 to SU{2) indicates 
the existence of novel unstable solutions. In Euclidean spacetime, it hints at a bar- 
rier between zero-action configurations in the trivial and non-trivial sectors (when 
spacetime is compactified to 5*4) - the /*. The Euchdean action would have a saddle 
point at the /* and it could significantly contribute to the path integral. When we 
consider static configurations, then the map suggests the existence of a solution with 
two directions of instability - the S*. Quarks and leptons are expected to have zero 
modes in the background of this configuration, which makes it a promising candidate 
for a fermionic soliton that allows a heavy fermion to decouple from the theory. Re- 
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stricting static configurations to be trivial in one spatial direction, the map suggests 
the existence of string solutions with three directions of instabihty - the W — string*. 
However, wc have not succeeded in constructing these solutions. The topological 
prescription suggests their existence and points to the region in configuration space 
where we can look for them. But it does not guarantee that the solutions exist. 

In 1 + 1 dimensional 0^ theory, we have reviewed the existence of approximate 
breathers - long-lived configurations that are localized in space and oscillate in time. 
Their existence hinges on two crucial properties of the theory: (1) there should be no 
massless particles in the theory; (2) the field potential should be non-linear in such 
a way that large amplitude oscillations would have a frequency lower than the lowest 
linear frequency (the lightest particle mass). When these two conditions are met, a 
configuration with fundamental frequency below the lowest linear frequency is stable 
against linear decay. However, it is not absolutely stable because the non-linearity 
means that higher harmonics of the fundamental frequency are present, and these 
couple to the linear modes. Nevertheless, such approximate breather configurations 
have lifetimes many orders of magnitude larger than all scales in the theory, which 
presents a challenge to the notion of naturalness in field theory. The Higgs-gauge 
sector of the electroweak theory has both properties (when the hypercharge gauge 
fields are ignored so that there is no unbroken subgroup and all three W bosons are 
massive) required for the existence of breathers. So we expect approximate breathers 
in the theory, and we are in the process of looking for them within a spherical ansatz. 
Their physical import could lie in the fact that during the electroweak phase tran- 
sition they would create out-of-equilibrium regions in space and thus contribute to 
electroweak baryogenesis. 

We have looked for stable solitons in the electroweak theory that carry the quan- 
tum numbers of a heavy fermion doublet and could allow it to decouple. We consider 
Higgs-gauge configurations within a spherical ansatz and compute fermion vacuum 
fiuctuation contributions to the energy, but ignore those of the bosonic fields. We find 
significant quantum corrections to the height of the sphaleron barrier. As we make 
the fermion heavier than the quantum-corrected sphaleron, we see the emergence of 
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a new barrier maintaining the exponential suppression of the fermion decay. For 
even larger values of the Yukawa coupling, however, the barrier disappears and the 
fermion's decay is unsupprcsscd. Wc do not sec any evidence for a fermionic soliton, 
and find that the fermion vacuum energy destabilizes would-be solitons. However, it 
is still possible, indeed likely, that the fermionic soliton that maintains anomaly can- 
cellation in the low energy theory and allows a heavy fermion to decouple lies outside 
the spherical ansatz. A most promising candidate is the -S"*, which is expected to exist 
due to the same non-trivial topology that gives rise to Witten's global anomaly. How- 
ever, in the absence of sufficient symmetry, the calculation of the one-loop effective 
energy required to analyze its stability becomes intractable. 

We have argued that it is hkely that classically unstable electroweak strings may 
be stabilized when quantum effects of quarks are considered, especially when several 
quarks are trapped along the string. Such stable multi-quark objects would have rich 
phenomenology associated with them, especially with regard to electroweak baryoge- 
nesis and dark energy. Such backgrounds generate long-range potentials for fermions, 
thereby presenting difficulties with our method to compute one-loop energies. So we 
perform the calculation for a similar, but simpler, problem of magnetic flux tubes 
in QED, as a stepping stone to the electroweak strings calculation. We find that 
when magnetic flux tubes are embedded in a no-net-flux conflguration (by including 
a spread-out region of return flux), the potential becomes short-range, and all sub- 
tleties and puzzles disappear. Moreover, as the return flux is made more diffuse, the 
energy approaches a well-deflned limit which corresponds to the energy of an isolated 
flux tube. This suggests that a similar embedding may be useful in the electroweak 
strings calculation. We also flnd that the quantum corrections to the energy in two 
and three spatial dimensions are very similar when we impose the same renormaliza- 
tion conditions in the two cases (which means a flnite renormalization in two spatial 
dimensions). This has lead us to currently investigate whether electroweak strings in 
two spatial dimensions can be stabilized by fermion quantum effects. This should be 
relevant for the stability in the physical case of three spatial dimensions. 

Undoubtedly, exploring physics beyond the Standard Model is intriguing and crit- 
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ical. But, as argued in this work, there are several unexplored objects, with wide- 
ranging significance, within the Standard Model itself. They address notions of nat- 
uralness and decoupling of fermions in chiral gauge theories. On a phenomenological 
front, they provide viable scenarios for electroweak baryogenesis and dark energy. We 
have explored some of these possibilities and much still remains to be discovered. 
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Appendix A 

Feynman Diagrams and 
Phaseshifts in the Spherical Ansatz 



A.l Results from Feynman Diagrams 

Here we list the results from the Feynman diagram calculations mentioned in Sec. 13.21 
In dimensional regularization {d = A — e) the counterterm coefficients, defined in 
eq. (jSHH), read 
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We have introduced the abbreviations 
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(A.2) 



where /i is the momentum scale introduced to maintain the canonical dimensions of 
the parameters when regularizing in fractional dimensions. 
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In eq. fj3.21|) E^^'"^^ denotes the contribution to the vacuum polarization energy 
from first and second order renormahzed Feynman diagrams. Its explicit expression 
reads 
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(A.3) 



with the Fourier transform of a field v^(a;) defined in the usual way as (p{q) = 
J d^xip{x)e^'^'^ . The third and fourth order counterterm contribution combined with 
the divergences in the third and fourth order Feynman diagrams is 
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where = $ — parameterizes the deviation of the Higgs field from its vev. 



A. 2 The Dirac Equation 

In this section, we describe how we obtain the bound state energies of the Dirac 
equation and the scattering phase shifts (and their Born series) in the presence of a 
background potential. These quantities are required to compute the vacuum polar- 
ization energy in eq. ()3.21|) . 

The fermion field obeys the time-independent Dirac equation 

Hd"^ = cu^ , (A.5) 

where 

Hd = -«7°f -9, + 7° K + V{^, W,)] , (A.6) 

and V is given in eq. ()3.7p . It is most convenient to use the chiral representation of 
the Dirac matrices, 



mf{s — ipTjXj) —iajdj 



The grand spin G is defined as the vector sum of isospin, spin and orbital angular 
momentum. It commutes with Ho as long as the fermion doublet is degenerate in 
mass and the background fields are in the spherical ansatz. We satisfy both conditions. 
For a given grand spin quantum number G (we suppress the grand spin projection 
label M throughout), the Dirac spinor has eight components and may be written 
in terms of generalized spherical harmonic functions 3^j,i(x) with j = G ± ^ and 
/ = j ± I as 

^G\X) =1 ' ' I , (A. 

«/l3^G+i G+1 + /23^G+i G + /sJ^G-i G + «/43^G-i G-1 
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where Qi and fi are radial functions and we have suppressed the grand spin labels on 
them. Note that in this chiral theory modes of different parity, e.g. gi and g2 mix. 
The spherical harmonics are two-component spinors in both spin and isospin space. 
The special case \E'o is defined only in terms of ^ and q and does not contain 
9^1 9a, /s) /4- 

The matrix elements of operators like TjXj between the spherical harmonics may 
be found in the literature ^U]- We use them to write the Dirac equation ()A.6|) as a 
a set of eight coupled first-order linear differential equations in the radial functions, 
for fixed G. From these equations we obtain the bound state solutions {\uj\ < rrif) 
in each grand spin channel using shooting algorithms. From Levinson's theorem we 
determine the number, iV^°"°*^, of bound states to shoot for, using phase shifts, 6g{uj), 
of the scattering state solutions of the Dirac equation: 

^bound ^ 1 ) _ ^ Sci-mj) - (5g(-oo)) . (A.9) 

To construct these scattering state solutions we re-write the Dirac equation as a set 
of second-order differential equations in the radial functions. Formally they read. 



hi2h2i-h^2{h22-^)K2{hii-u) ^g = 0, (A.IO) 



where 



= myG+^,G+l + ^23^G+i G + ^33^G-i G + ^^43^G-i G-1 



denotes the upper two-component spinor in ^g- Iii fhe chiral representation of the 
Dirac matrices, we require s^(r) +p^(r) > so that hi2 is invertible. As mentioned 
in Sec. 13.4.21 this is a restriction on our variational ansatze. Using the known matrix 
elements for the spin-isospin operators like TjXi, we then project eq. (jA.lOj) onto grand 
spin channels and obtain the desired second order differential equations. These are 
rather lengthy and are listed a the end of this Appendix. They may be written in the 
form, 

^|DG(r) + ArG(r)|- + MG(r)| ^?,(r) = 0, (A.ll) 
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with 





Og 



diag ((G + 1)(G + 2), G{G + 1), G{G + 1), (G - 1)G) 



(A.12) 



and k"^ = uj'^ — m'j. The matrices Ndr) and MG'(r) are given in terms of the functions 
s(r), . . . ,7(r) that specify the static background fields in the spherical ansatz, as in 
eq. (j3.26|) . Their elements can be read off from the equations displayed at the end of 
this Appendix. As r — » oo, Ndr) and Mdr) and the differential equations 
decouple, as long as the potential goes to zero sufficiently fast. 

We have a four-channel scattering problem. We express the four wavefunctions 
and four boundary conditions in matrix form, Gij{r) = g\^\r), where the linearly 
independent boundary conditions are labeled by j = 1,2,3,4. We then write Q{r) 
as a multiplicative modification of the matrix solution to the free differential equa- 
tions, Q{r) = F{r)-H{kr), where H{x) = diag{hQ\_i{x) , hQ\x), hQ\x), h^GL^{x)) with 
h^l\x) denoting spherical Hankel functions of the first kind such that Dcir) ■ H{kr) = 
0. (The 4x4 matrices F and H depend on the grand spin quantum number G. For 
convenience we omit that label from now on.) Imposing the boundary conditions 
F{r — s> oo) = and F'{r — * oo) = 0, it is clear that the i*'^ row of Q describes an 
outgoing spherical wave in the i^^ channel. Similarly, Q* describes incoming spherical 
waves. The scattering wavefunction can be written as 



g*{r) + g{r)S{k), 



(A. 13) 



and requiring this to be regular at the origin gives the scattering matrix 



S{k) = \im H-\kr)F~\r)F*{r)H*{kr) . 



(A. 14) 
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We are interested in the sum of the eigenphase shifts in a given grand spin channel, 
S{k) = —Tiln S{k) = — limTrln(F^i(r)F*(r)) . (A.15) 



An efficient way to avoid any ambiguities in additive contributions of multiples of vr 
in 6{k) is to define 

5{k,r) = —TT\n(^F-\r)F*{r)) , (A.16) 



2i 

with 6{k) = 6{k,0). We then integrate 



-^^Tr F'F'M (A. 17) 



dr 

along with F{k, r) from infinity to with the boundary condition limj,^oo ^{k, r) = 
to obtain 6{k) as a smooth function of k. The differential equation for the matrix 
F{k,r), 

= F" + -F' + 2F'L' + \ [F, O] + N{F' + FL') + MF , (A. 18) 



is obtained from 



D{r) + N{r)^ + M{r)\g{r) 



H-\kr)=0, (A.19) 



where L{kr) = In H{kr) and primes denote derivatives with respect to the radial 
coordinate. The components of L'{kr) can be expressed as simple rational functions, 
which avoids any numerical instability that would be caused by the oscillating Hankel 
functions. 

To construct the Born series for S{k), we introduce F^^\k,r) where n labels the 
order in the background fields in an expansion around the vacuum configuration with 
f^^^(r) = in eq. 1)3.311) . We obtain the corresponding differential equations 



= + Vi)' + 2F«'L' + l[F«,0] + Ar«L' + M(i), 

= F(2)%^F(2)' + 2F(2)'L' + l[F(2),0] + Ar«(F«' + F«L') 

+Ar(2)L' + M(i)f(i) + M(2) , (A.20) 
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where the matrices A^*^*^ and M^^^ are obtained from and M by expanding to 
order i in the deviation of the background fields from the above described vacuum 
configuration. We integrate these differential equations with the boundary conditions 
F^^\k, oo) = and oo) = and obtain 



We eliminate the quadratic divergence from the vacuum polarization energy by sub- 
tracting these from 6{k) and adding them back in as renormalized first and second 
order Feynman diagrams. 

There still remains the logarithmic divergence whose elimination would require 
third and fourth order Born subtractions. These become considerably more compli- 
cated, so instead we use the limiting function approach as described in [21]. The idea 
is to subtract only the local contributions to the third and fourth Born approximants 
to the phase shift by identifying them with the divergent contributions to the third 
and fourth order Feynman diagrams. To this end we formally manipulate these di- 
vergent Feynman diagrams. To extract the local contributions we set the external 
momenta to zero and then integrate over the energy and the two spatial angles of the 
loop momenta, k^, such that a (regularized) integral over k = \k\ is left. We write its 
integrand in the form as in eq. ()3.2H) . 



S'-^\k) 
5^^\k) 



Str (f(^)(A;,0)) , 




(A.21) 




where 





(A.22) 
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in the Wq = gauge and where = $ — f denotes the deviation of the Higgs field 
from its vev. 

Thus we numerically determine the bound state energies, the phase shifts, their 
Born series, and the limiting function. These are all ingredients in the expression for 
the fermion vacuum polarization energy in eq. (j3.2H) . 
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A. 2.1 The Second Order Equations 

Here we list the four coupled second order equations in the four radial functions 
specifying the left-handed spinor for a fixed grand spin G, obtained from the Dirac 
equation. The matrix elements of Mq, Nq in eg. lA. Ill can be read off from these, and 
then the phaseshifts computed, as described above. The radial functions a,7,ai,S 
and T] specify static background Higgs-gauge fields in the spherical ansatz, as in 
eqs. KM EM 







9i +9i 
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,r2 1 di:' 1 + G . , 

+ . (2 sm [j]) - 7) 



(l + 2G')r 



r 2S2 dr (l + 2G)r 



--flir - (1 + G)a - rr( ^{\^G) sin(2?7) 
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r{l + 2G) 
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Appendix B 



The QED Flux Tube Phaseshifts 



We explain how fermion phaseshifts and their Born series are computed in QED in 
the background of magnetic flux tubes. These are ingredients in the fermion quantum 
corrections to the energy of flux tubes, as discussed in Sec. 14.41 

B.l The Dirac Equation 

Consider QED in 2+1 dimensions with a four-component spinor. In the radial gauge, 
vortex conflgurations of gauge flelds give flux F tubes: 

Ao = 0, 

A ^ ^/(p)^. (B.l) 

For small p, /(p) oc p^, otherwise B is singular. For p large compared to the width of 
the flux tube, /(p) approaches 1. In the case of no-net-flux embedding by considering 
a distant region of return flux, /(p) approaches for p beyond the region of return 
flux. The energy uj eigenspinor, obeys the time-independent Dirac equation 

Hd^ = u;^ . (B.2) 
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The Dirac Hamiltonian is 



Hd = -i 



+ 7° (m + F) , 



(B.3) 



where the potential is 



V = - — /(p) (^-7^ sin + 7^ cos 



(B.4) 



Note that in the case of a flux tube with no return flux, the potential falls slowly hke 
1/p at large p, and conflicts with standard assumptions of scattering theory. However, 
in the no-net-flux embedding, the potential is short-range. 

Since the vortex configuration is cylindrically symmetric, the total angular mo- 
mentum Jz commutes with the above Hamiltonian. The subscript z is meaningless in 
2-1-1 dimensions, but we include it because we are also interested in 3-1-1 dimensions 
in which case the z projection of the total angular momentum is a good quantum 
number. For a given quantum number (say M, which is a half integer), the Dirac 
spinor has the form 



where we have suppressed the label on the radial functions /j, gi. The kets |+), |— ) 
are spin up and spin down eigenkets of Sz- In the chiral representation of the Dirac 
matrices (in which 75 is diagonal), the upper and lower components of the above 
spinor denote left and right handed projections respectively. 

We use the first order Dirac equation to eliminate the left handed spinor and 
obtain second order differential equations for the radial functions of the right handed 
spinor: 




(B.5) 












(B.6) 
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where k"^ = u'^ — m?, the primes denote derivatives with respect to p, f is the profile 
function specifying the vortex configuration, and = ff • 

B.2 Phaseshifts 
B.2.1 Zero Flux Case 

First consider the no-net-flux embedding in which /(p oo) = 0. The asymptotic 
equations for large p have outgoing wave solutions that are Hankel functions of the 
first kind. We write the solutions in the presence of the background as a multiplicative 
modification of the free solutions: 

g,{p) = e^^^^^^Hll^^ikp),g,ip) = e^'^^('')/f«^,|(M • (B.7) 

Then, plugging in this form into the differential equations, we get the equations for 
A: 

= tl3'; + 2iL',j_.J[ + - {^[f - ^ [2 (M - i) / + p/' + :Ff] , 

= i^'l + 2^/5^L;,^. + 1/5^ - {P',f - ^ [2 (M + i) / - p/' + Tf] , (B.8) 

where Ly{kp) = In H\^\{kp). The boundary conditions at infinity are that = = 0. 
Integrating the (3i from oo to 0, the total phaseshift for = M is (summing over the 
two scattering channels and over the two signs of the energy) 

5M(A;) = -2Re(/3i(0)+/32(0)) . (B.9) 

(See j2ni for a derivation of the relationship between phaseshifts and (3.) 

To obtain a Born series expansion of the phaseshifts, we introduce where j 
labels the order in an expansion around A = 0. The equations for the first two terms. 
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/3f\/3f\are 










))" + 2^L'^_i (tf V + -SP?)' - ^ [2 (M - i) / + pf] , 



P P 



P 



P' 



P' 



with vanishing /Sp-* and {/S!-^')' as p — > oo. Then, the first two terms in the Born series 
expansion of the phaseshifts are obtained by integrating the to the origin: 



(B.IO) 



(B.ll) 



B.2.2 Non-zero Flux case 

In the case of non-zero flux, the profile function /(p) goes to 1 at large p. This gives 
a long-range scattering potential. As before, expressing the radial functions gi{p) as 
multiplicative modifications of the asymptotic solutions (Hankel functions with their 
order shifted by the flux). 



(B.12) 



Then, plugging in this form into the differential equations, we get the equations for 
A: 







+- 



T r 



2(M-i)(l-/) + +^(l-/2)-/V 



= 



i/3^' + 2iL'^^.^^/3^ + -/3^-(/3^)2 



2(M+i)(l-/) + +^(l-/2)+/V 



(B.13) 
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where, as before, L^{kp) = \nH\^\{kp), and at infinity, = = 0. Integrating tlie 
Pi from oo to 0, the total phaseshift for = M is (summing over the two scattering 
channels and over the two signs of the energy) 

OM, ideal (fc) = -2Re(/5i(0)+/?2(0)) . (B.14) 

However, as denoted by the subscript, the above phaseshift is relative to an ideal 
(0 width) flux tube, because the asymptotic equations do not correspond to the free 
equations. So, as discussed in further detail in Sec. 14.4.31 the phaseshift of the ideal 
vortex relative to the totally-trivial (A^ = 0) configuration must be added to the 
above result: 

SMMeik) = 7r(|M-i| + |M+i|-|M-i+^|-|M+i+^|) . (B.15) 

Although we could generate the Born series expansion of the phaseshift using an 
expansion of the /3j as before, there is the complication that the asymptotic solutions 
are not free solutions because of the long-range potential. This can be dealt with, of 
course. However, we use the fact that the Born series expansion is an expansion in 
powers of the flux, JF, and we compute the terms using numerical derivatives of the 
phaseshift with respect to the flux evaluated at JF = 0: 



2 rfj^2 



(B.16) 

:f=o 



Thus, we compute the phaseshifts and the first two terms in their Born series, 
which allows us to exactly compute the renormalized one-loop energies of magnetic 
flux tubes in QED in Sec. EM 
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